Controller and control method for internal combustion engine

ABSTRACT

To provide a controller and a control method for an internal combustion engine capable of reducing the calculation error of recirculation exhaust gas amount due to changes with time of the internal combustion engines, and humidity change of intake air, and also capable of reducing the calculation error of recirculation exhaust gas amount at transient operation. The controller and the control method for the internal combustion engine calculates humidity detecting EGR rate based on intake-air humidity and manifold humidity, calculates humidity detecting opening area which realizes humidity detecting recirculation flow rate calculated based on humidity detecting EGR rate, calculates learned opening area corresponding to present opening degree of EGR valve using learning value of opening area calculated based on humidity detecting opening area, and calculates flow rate of recirculation exhaust gas for control based on learned opening area.

INCORPORATION BY REFERENCE

The disclosure of Japanese Patent Application No. 2016-67075 filed on Mar. 30, 2016 including its specification, claims and drawings, is incorporated herein by reference in its entirety.

BACKGROUND

The present invention relates to a controller and a control method for an internal combustion engine that is provided with an intake path and an exhaust path, a throttle valve for opening and closing the intake path, an EGR path for recirculating exhaust gas from the exhaust path to an intake manifold that is part of the intake path and is at the downstream side of the throttle valve, and an EGR valve for opening and closing the EGR path.

In order to appropriately control an internal combustion engine, it is important to accurately calculate the amount of air taken into a cylinder and to accurately control the fuel supply amount and the ignition timing in accordance with the cylinder intake air amount. Ignition timing control needs to change an ignition timing to the ignition timing (MBT: Minimum Spark Advance for Best Torque) at which the output torque becomes maximal or the like, in accordance with not only the rotation speed of the internal combustion engine and the cylinder intake air amount but also other factors such as the coolant temperature of the internal combustion engine, whether or not a knock has occurred, the fuel property, and the EGR (Exhaust Gas Recirculation) rate.

Meanwhile, with regard to the EGR, there are two methods, i.e., a method (referred to as an external EGR, hereinafter) in which an EGR valve is provided in an EGR path for recirculating exhaust gas from the exhaust path to the intake manifold and the amount of exhaust gas to be recirculated to the intake manifold is controlled based on the opening degree of the EGR valve and a method (referred to as an internal EGR, hereinafter) in which a variable valve timing mechanism, which makes the opening/closing timings of one of or both of an intake valve and an exhaust valve variable, is provided and an overlap period, during which the intake valve and the exhaust valve are concurrently opened, is changed so that the amount of exhaust gas remaining in the cylinder is controlled. In recent years, in order to reduce the fuel cost and raise the output, the number of internal combustion engines provided with both the external EGR mechanism and the internal EGR mechanism has been increasing. In the present invention, an EGR and an EGR rate, when simply described in this manner, denote an external EGR and an external EGR rate, respectively.

In recent years, an internal combustion engine has been controlled by utilizing, as an index, the output torque of the internal combustion engine. Because the thermal efficiency changes depending on the cylinder intake air amount and the EGR rate, it is required to estimate the output torque based on the cylinder intake air amount and the EGR rate. Accordingly, either in order to control the ignition timing or in order to estimate the output torque, it is required to accurately estimate the EGR rate.

As a technology for estimating an EGR rate, for example, the technologies disclosed in Japanese Patent No. 5642222 and Japanese Examined Patent Publication No. S58-55345 have already been known. In the technology disclosed in Japanese Patent No. 5642222, by use of the intake air amount detected by an air flow sensor, the cylinder flow rate, which is the amount of air that flows into the cylinder and is calculated based on the pressure in the intake manifold or the like, and the opening degree of the EGR valve, the variation in the characteristics of the EGR valve and the change with time thereof are learned so that the recirculation exhaust gas amount is estimated.

The technology disclosed in Japanese Examined Patent Publication No. S58-55345 is to perform feedback control of the opening degree of an EGR valve so that the CO2 (carbon dioxide) concentration detected by a CO2 concentration sensor provided in the intake manifold approaches a target value. Japanese Examined Patent Publication No. S58-55345 also discloses a configuration in which instead of a CO2 concentration sensor, a humidity sensor is provided.

SUMMARY

In the technology disclosed in Japanese Patent No. 5642222, it is not required to add a dedicated sensor for estimating a recirculation exhaust gas amount and hence the cost hike is not caused by the increase in the number of components; however, because the recirculation exhaust gas amount is indirectly estimated, there has been a problem that an estimation error is caused by the individual differences, the changes with time, and the like in the respective characteristics of internal combustion engines and various kinds of sensors.

In the technology disclosed in Japanese Examined Patent Publication No. S58-55345, the recirculation exhaust gas amount is feedback-controlled based on the humidity in the intake manifold; however, because the effect of the humidity of intake air that is newly taken into the intake manifold from the atmosphere, i.e., the effect of the humidity of the atmospheric air, which is provided to the humidity in the intake manifold, is not taken into consideration, there has been a problem that an error in controlling the recirculation exhaust gas amount is caused. Because the humidity of the atmospheric air largely changes depending on the district, the season, the weather, and the like, the error in controlling the recirculation exhaust gas amount exceeds the upper limit that can be neglected.

A general humidity sensor has a response delay of time constant of second order. Therefore, there has been a problem that a response delay caused in the recirculation exhaust gas amount calculated based on the output signal of the humidity sensor, and the calculation error of the recirculation exhaust gas amount at the time of transient operation becomes large.

Thus, there has been desired a controller and a control method for an internal combustion engine capable of reducing the calculation error of the recirculation exhaust gas amount due to the individual differences and the changes with time of the internal combustion engines, and the humidity change of intake air (atmospheric air), and also capable of reducing the calculation error of the recirculation exhaust gas amount at the time of transient operation.

According to a first aspect of the present invention, a controller for an internal combustion engine that is provided with an intake path and an exhaust path, a throttle valve for opening and closing the intake path, an EGR path for recirculating exhaust gas from the exhaust path to an intake manifold that is part of the intake path and is at the downstream side of the throttle valve, and an EGR valve for opening and closing the EGR path, the internal combustion engine controller includes:

a driving-condition detector that detects a manifold pressure, which is a pressure of gas in the intake manifold, a manifold temperature, which is a temperature of gas in the intake manifold, a manifold humidity, which is a humidity of gas in the intake manifold, an intake-air pressure, which is a pressure of intake air to be taken into the intake path, an intake-air temperature, which is a temperature of the intake air, an intake-air humidity, which is a humidity of the intake air, an intake air flow rate, which is a flow rate of the intake air, and an opening degree of the EGR valve;

a humidity detecting EGR rate calculator that calculates a humidity detecting EGR rate, which is a ratio of a recirculation exhaust gas, which is the exhaust gas recirculated into the intake manifold, to the intake air, based on the intake-air temperature, the intake-air humidity, the intake-air pressure, the manifold temperature, the manifold humidity, and the manifold pressure;

an opening area learning value calculator that calculates a humidity detecting recirculation flow rate which is a flow rate of the recirculation exhaust gas based on the humidity detecting EGR rate and the intake air flow rate, calculates a humidity detecting opening area, which is an opening area of the EGR valve which realizes the humidity detecting recirculation flow rate, and calculates a learning value of the opening area of the EGR valve based on the humidity detecting opening area; and

a recirculation exhaust gas calculator for control that calculates a learned opening area of the EGR valve corresponding to the present opening degree of the EGR valve using the learning value of the opening area, and calculates a flow rate of the recirculation exhaust gas for control used for controlling the internal combustion engine based on the learned opening area.

According to a second aspect of the present invention, a control method for an internal combustion engine is a control method for an internal combustion engine that is provided with an intake path and an exhaust path, a throttle valve for opening and closing the intake path, an EGR path for recirculating exhaust gas from the exhaust path to an intake manifold that is part of the intake path and is at the downstream side of the throttle valve, and an EGR valve for opening and closing the EGR path, the control method includes:

a driving-condition detecting that detects a manifold pressure, which is a pressure of gas in the intake manifold, a manifold temperature, which is a temperature of gas in the intake manifold, a manifold humidity, which is a humidity of gas in the intake manifold, an intake-air pressure, which is a pressure of intake air to be taken into the intake path, an intake-air temperature, which is a temperature of the intake air, an intake-air humidity, which is a humidity of the intake air, an intake air flow rate, which is a flow rate of the intake air, and an opening degree of the EGR valve;

a humidity detecting EGR rate calculating that calculates a humidity detecting EGR rate, which is a ratio of a recirculation exhaust gas, which is the exhaust gas recirculated into the intake manifold, to the intake air, based on the intake-air temperature, the intake-air humidity, the intake-air pressure, the manifold temperature, the manifold humidity, and the manifold pressure;

an opening area learning value calculating that calculates a humidity detecting recirculation flow rate which is a flow rate of the recirculation exhaust gas based on the humidity detecting EGR rate and the intake air flow rate, calculating a humidity detecting opening area, which is an opening area of the EGR valve which realizes the humidity detecting recirculation flow rate, and calculates a learning value of the opening area of the EGR valve based on the humidity detecting opening area; and

a recirculation exhaust gas calculating for control that calculates a learned opening area of the EGR valve corresponding to the present opening degree of the EGR valve using the learning value of the opening area, and calculates a flow rate of the recirculation exhaust gas for control used for controlling the internal combustion engine based on the learned opening area.

According to the controller and the control method for the internal combustion engine, regardless of the individual differences and the changes with time of the internal combustion engines, and the humidity change of intake air, the EGR rate (the humidity detecting EGR rate) can be detected with sufficient accuracy based on the intake-air humidity and the manifold humidity. Then, based on the humidity detecting EGR rate, the learning value of the opening area of the EGR valve is calculated. Therefore, even if it is a case where the flow characteristic of the EGR valve changes by deposit, such as soot, and a case where an EGR valve stops operating by aging degradation, the flow characteristic of the EGR valve can be learned with sufficient accuracy. Then, based on the learned opening area corresponding to the present opening degree of the EGR valve, the flow rate of the recirculation exhaust gas for control can be calculated with good responsiveness. Therefore, while reducing the individual differences and the changes with time of the internal combustion engines, and the humidity change of intake air, the calculation error of the recirculation exhaust gas amount at the time of transient operation can be reduced, and the control accuracy of the internal combustion engine can be improved.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic configuration diagram of an internal combustion engine and a controller according to Embodiment 1 of the present invention;

FIG. 2 is a block diagram of a controller according to Embodiment 1 of the present invention;

FIG. 3 is a hardware configuration diagram of the controller according to Embodiment 1 of the present invention;

FIG. 4 is a chart representing the state of partial pressures of gasses in the intake manifold according to Embodiment 1 of the present invention; and

FIG. 5 is a flowchart representing the processing by the controller according to Embodiment 1 of the present invention.

DETAILED DESCRIPTION OF THE EMBODIMENTS 1. Embodiment 1

A controller 50 for an internal combustion engine 1 (hereinafter, referred to simply as the controller 50) according to Embodiment 1 will be explained with reference to the drawings. FIG. 1 is a schematic configuration diagram of the internal combustion engine 1 and the controller 50; FIG. 2 is a block diagram of the controller 50 according to Embodiment 1. The internal combustion engine 1 and the controller 50 are mounted in a vehicle; the internal combustion engine 1 functions as a driving-force source for the vehicle (wheels).

1-1. The Configuration of the Internal Combustion Engine 1

The configuration of the internal combustion engine 1 will be explained. The internal combustion engine 1 has a cylinder 25 in which a fuel-air mixture is combusted. The internal combustion engine 1 has an intake path 23 for supplying air to the cylinder 25 and an exhaust path 17 for discharging exhaust gas combusted in the cylinder 25. The internal combustion engine 1 has a throttle valve 6 for opening and closing the intake path 23. The throttle valve 6 is an electronically controlled throttle valve that is opening/closing-driven by an electric motor controlled by the controller 50. The throttle valve 6 is provided with a throttle opening degree sensor 7 that generates an electric signal according to a throttle opening degree of the throttle valve 6.

An air cleaner 24 for purifying air taken into the intake path 23 is provided at the most upstream portion of the intake path 23. In the intake path 23 at the upstream side of the throttle valve 6, there are provided an air flow sensor 3 that outputs an electric signal according to the flow rate of intake air, which is air to be taken from the atmosphere into the intake path 23, an intake-air temperature sensor 4 that outputs an electric signal according to an intake-air temperature Ta, which is the temperature of intake air, and an intake-air humidity sensor 5 that outputs an electric signal according to an intake-air humidity Hra, which is the humidity of intake air. The pressure in the intake path 23 at the upstream side of the throttle valve 6 can be regarded as equal to the atmospheric pressure. An intake-air pressure sensor 2 that outputs an electric signal according to an intake-air pressure Pa, which is the pressure of intake air (the atmospheric air, in this example), is provided outside the intake path 23 (for example, inside the controller 50).

Each of the intake-air temperature sensor 4 and the intake-air humidity sensor 5 may be either integrated with or separated from the air flow sensor 3. Alternatively, each of the intake-air temperature sensor 4 and the intake-air humidity sensor 5 may be provided outside the intake path 23, as is the case with the intake-air pressure sensor 2; the intake-air pressure sensor 2 may be provided at a place the same as the place where the intake-air temperature sensor 4 and the intake-air humidity sensor 5 are provided. In any case, the intake-air pressure sensor 2, the intake-air temperature sensor 4, and the intake-air humidity sensor 5 are provided at a place where there exists intake air to be taken into the intake path 23 and the pressure of the intake air is substantially the same.

The portion, at the downstream side of the throttle valve 6, of the intake path 23 is an intake manifold 12. The upstream portion of the intake manifold 12 is a surge tank 11 for suppressing an intake-air ripple. The internal combustion engine 1 has an EGR path 21 for recirculating exhaust gas from the exhaust path 17 to the intake manifold 12 (the surge tank 11, in this example) and an EGR valve 22 for opening and closing the EGR path 21. The EGR valve 22 is an electronically controlled EGR valve that is opening/closing-driven by an electric actuator such as an electric motor to be controlled by the controller 50. The EGR valve 22 is provided with an EGR opening degree sensor 27 that outputs an electric signal according to the opening degree Oe of the EGR valve 22. Exhaust gas (referred to as recirculation exhaust gas, hereinafter) that is recirculated to the surge tank 11 and intake air that is taken into the surge tank 11 are mixed and homogenized in the surge tank 11. “EGR” is an acronym for Exhaust Gas Recirculation.

In the intake manifold 12, there are provided a manifold pressure sensor 8 that outputs an electric signal according to a manifold pressure Pb, which is the pressure of gas in the intake manifold 12, a manifold temperature sensor 9 that outputs an electric signal according to a manifold temperature Tb, which is the temperature of gas in the intake manifold 12, and a manifold humidity sensor 10 that outputs an electric signal according to a manifold humidity Hrb, which is the humidity of gas in the intake manifold 12. Each of the manifold temperature sensor 9 and the manifold humidity sensor 10 may be either integrated with or separated from the manifold pressure sensor 8. The manifold temperature sensor 9 and the manifold humidity sensor 10 are provided at the downstream side of the connection point between the intake manifold 12 and the EGR path 21 so that the temperature and the humidity of a gas obtained by sufficiently mixing intake air with recirculation exhaust gas can be detected.

An injector 13 for injecting a fuel is provided at the downstream side part of the intake manifold 12. The injector 13 may be provided so as to inject a fuel directly into the cylinder 25.

An ignition plug for igniting a fuel-air mixture and an ignition coil 16 for supplying ignition energy to the ignition plug are provided on the top of the cylinder 25. On the top of the cylinder 25, there are provided an intake valve 14 for adjusting the amount of intake air to be taken from the intake path 23 into the cylinder 25 and an exhaust valve 15 for adjusting the amount of exhaust gas to be exhausted from the cylinder to the exhaust path 17. The intake valve 14 is provided with an intake variable valve timing mechanism that makes the opening and closing timing thereof variable. The exhaust valve 15 is provided with an exhaust-gas variable valve timing mechanism that makes the opening/closing timing therefor variable. Each of the variable valve timing mechanisms 14 and 15 has an electric actuator. On the crankshaft of the internal combustion engine 1, there is provided a crank angle sensor 20 for generating an electric signal according to the rotation angle thereof.

In the exhaust path 17, there is provided an air-fuel ratio sensor 18 for detecting an air-fuel (Air/Fuel) ratio AF, which is the ratio of air to fuel in exhaust gas. A catalyst 19 for purifying exhaust gas is also provided in the exhaust path 17. As the catalyst 19, a three-way catalyst, the purification performance of which becomes higher in the vicinity of the theoretical air-fuel ratio AF0, is utilized.

1-2. The Configuration of the Controller 50

Next, the controller 50 will be explained. The controller 50 is the one whose control subject is the internal combustion engine 1. As represented in FIG. 2, the controller 50 is provided with control units such as a driving-condition detection unit 51, a humidity detecting EGR rate calculation unit 52, an opening area learning value calculation unit 53, a recirculation exhaust gas calculation unit for control 54, and a recirculation amount utilization control unit 55. The respective control units 51 through 55 and the like of the controller 50 are realized by processing circuits included in the controller 50. Specifically, as illustrated in FIG. 3, the controller 50 includes, as processing circuits, a computing processing unit (computer) 90 such as a CPU (Central Processing Unit), storage apparatuses 91 that exchange data with the computing processing unit 90, an input circuit 92 that inputs external signals to the computing processing unit 90, an output circuit 93 that outputs signals from the computing processing unit 90 to the outside, a communication circuit 94 in which the computing processing unit 90 performs data communication with external apparatuses, and the like.

As the storage apparatuses 91, there are provided a RAM (Random Access Memory) that can read data and write data from the computing processing unit 90, a ROM (Read Only Memory) that can read data from the computing processing unit 90, and the like. The input circuit 92 is connected with various kinds of sensors and switches and is provided with an A/D converter and the like for inputting output signals from the sensors and the switches to the computing processing unit 90. The output circuit 93 is connected with electric loads and is provided with a driving circuit and the like for outputting a control signal from the computing processing unit 90. The communication circuit 94 is connected with other vehicle electronic apparatuses such as an air conditioner controller 80, transmission controller 81, and the like through a communication wire and performs cable communication based on a communication protocol such as the CAN (Controller Area Network).

In addition, the computing processing unit 90 runs software items (programs) stored in the storage apparatus 91 such as a ROM and collaborates with other hardware devices in the controller 50, such as the storage apparatus 91, the input circuit 92, the output circuit 93, and the communication circuit 94, so that the respective functions of the control units 51 through 55 included in the controller 50 are realized. Setting data items such as characteristic data and determination values to be utilized in the control units 51 through 55 are stored, as part of software items (programs), in the storage apparatus 91 such as a ROM.

In Embodiment 1, the input circuit 92 is connected with the intake-air pressure sensor 2, the air flow sensor 3, the intake-air temperature sensor 4, the intake-air humidity sensor 5, the throttle opening degree sensor 7, the manifold pressure sensor 8, the manifold temperature sensor 9, the manifold humidity sensor 10, the air-fuel ratio sensor 18, the crank angle sensor 20, an accelerator position sensor 26, the EGR opening degree sensor 27, and the like. The output circuit 93 is connected with the throttle valve 6 (electric motor), the injector 13, the intake variable valve timing mechanism 14, the exhaust variable valve timing mechanism 15, the ignition coil 16, the EGR valve 22 (the electric actuator), and the like. The controller 50 is connected with various kinds of unillustrated sensors, switches, actuators, and the like.

As basic control, the controller 50 calculates a fuel injection amount, an ignition timing, and the like, based on inputted output signals and the like from the various kinds of sensors, and then performs driving control of the injector 13, the ignition coil 16, and the like. Based on the output signal of the accelerator position sensor 26 and the like, the controller 50 calculates the output torque of the internal combustion engine 1, demanded by the driver, and then controls the throttle valve 6 and the like so that an intake air amount for realizing the demanded output torque is obtained. In this case, the output torque of the internal combustion engine 1 which is described below and which is calculated based on the recirculation exhaust gas flow Qes for control may be taken into consideration. Specifically, the controller 50 calculates a target throttle opening degree and then performs driving control of the electric motor for the throttle valve 6 so that the throttle opening degree, detected based on the output signal of the throttle opening degree sensor 7, approaches the target throttle opening degree.

<The Driving-Condition Detection Unit 51>

The driving-condition detection unit 51 detects the driving conditions of the internal combustion engine 1 and the vehicle. The driving-condition detection unit 51 detects various kinds of driving conditions, based on, for example, the output signals of various kinds of sensors. The driving-condition detection unit 51 detects a manifold pressure Pb, a manifold temperature Tb, and a manifold humidity Hrb. In Embodiment 1, the driving-condition detection unit 51 detects the manifold pressure Pb, based on the output signal of the manifold pressure sensor 8. The driving-condition detection unit 51 detects the manifold temperature Tb, based on the output signal of the manifold temperature sensor 9. The driving-condition detection unit 51 detects the manifold humidity Hrb, based on the output signal of the manifold humidity sensor 10.

The driving-condition detection unit 51 detects an intake-air pressure Pa, an intake-air temperature Ta, and an intake-air humidity Hra. In Embodiment 1, the driving-condition detection unit 51 detects the intake-air pressure Pa, based on the output signal of the intake-air pressure sensor 2. The driving-condition detection unit 51 detects the intake-air temperature Ta, based on the output signal of the intake-air temperature sensor 4. The driving-condition detection unit 51 detects the intake-air humidity Hra, based on the output signal of the intake-air humidity sensor 5.

In Embodiment 1, as each of the intake-air humidity sensor 5 and the manifold humidity sensor 10, a sensor of the type that detects a relative humidity is utilized; for example, an electric-resistor sensor that detects a relative humidity based on the electric resistance value of a moisture-sensitive material, an electrostatic-capacitance sensor that detects a relative humidity based on the electrostatic capacitance of a sensor element, or the like is utilized. Accordingly, the driving-condition detection unit 51 detects a relative humidity, as each of the manifold humidity Hrb and the intake-air humidity Hra. The humidity sensors 5 and 10 have a response delay of the time constant for about several seconds.

The driving-condition detection unit 51 detects a throttle opening angle, based on the output signal of the throttle position sensor 7, and detects the opening degree Oe of the EGR valve 22, based on the output signal of the EGR opening degree sensor 27. The driving-condition detection unit 51 detects an air-fuel ratio AF of the exhaust gas, based on the output signal of the air-fuel ratio sensor 18, detects a rotational speed Ne of the internal combustion engine 1, based on the output signal of the crank angle sensor 20, and detects an accelerator opening degree, based on the output signal of the accelerator position sensor 26.

The driving-condition detection unit 51 detects the intake air flow rates Qa, based on the output signal of the air flow sensors 3. The driving-condition detection unit 51 calculates an intake air amount QA [g/stroke] taken into the intake path 23 (intake manifold 12) in one stroke period (for example, the interval of BTDC5degCA) based on the intake air flow rate Qa [g/s], as shown in an equation (1); and applies first-order-lag filtering processing, which simulates a delay in the intake manifold 12 (surge tank), to the intake air amount QA, so as to calculate a cylinder intake air amount QAc [g/stroke] taken into the cylinder 25 in one stroke period. The driving-condition detection unit 51 calculates the intake air amount QA by multiplying one stroke period AT to the intake air flow rate Qa, for example.

QAc(n)=KCCA·QAc(n−1)+(1−KCCA)·QA(n)

QA(n)=Qa(n)·ΔT(n)  (1)

Where KCCA is a preliminarily set filter gain. Where (n) denotes the value in the present calculation cycle, and (n−1) denotes the value in the immediately previous calculation cycle.

The driving-condition detection unit 51 calculates a charging efficiency Ec of intake air by dividing the cylinder intake air amount QAc by a value obtained by multiplying the density ρ0 of air under the standard atmospheric condition to the cylinder volume Vc. The charging efficiency Ec is the ratio of the cylinder intake air amount QAc to the mass (ρ0×Vc) of air under the standard atmospheric condition, with which the cylinder volume Vc is filled. The standard atmospheric condition denotes the state of 1 atm and 25° C.

$\begin{matrix} {{Ec} = \frac{QAc}{\rho \; {0 \cdot {Vc}}}} & (2) \end{matrix}$

The driving-condition detection unit 51 detects a temperature Tex of the exhaust gas at the exhaust path 17 side of the EGR valve 22. In Embodiment 1, the driving-condition detection unit 51 calculates the temperature Tex of the exhaust gas corresponding to the present rotational speed Ne of the internal combustion engine 1 and the present charging efficiency Ec, by use of an exhaust gas temperature characteristic data in which the relationship among the rotational speed Ne of the internal combustion engine 1, the charging efficiency Ec, and the temperature Tex of the exhaust gas. A temperature sensor may be provided in the exhaust path 17, and the driving-condition detection unit 51 may detect the temperature Tex of the exhaust gas based on the output signal of the temperature sensor. A data map, a data table, a polynomial, an equation, and the like are used for each characteristic data.

The driving-condition detection unit 51 detects a pressure Pex of the exhaust gas at the exhaust path 17 side of the EGR valve 22. In Embodiment 1, the driving-condition detection unit 51 calculates the pressure Pex of the exhaust gas corresponding to the present rotational speed Ne of the internal combustion engine 1 and the present charging efficiency Ec, by use of an exhaust gas pressure characteristic data in which the relationship among the rotational speed Ne of the internal combustion engine 1, the charging efficiency Ec, and the pressure Pex of exhaust gas.

<The Humidity Detecting EGR Rate Calculation Unit 52>

The humidity detecting EGR rate calculation unit 52 calculates a humidity detecting EGR rate Regr, which is a ratio of the recirculated exhaust gas, which is the exhaust gas recirculated into the intake manifold 12, to the intake air, based on the intake-air temperature Ta, the intake-air humidity Hra, the intake-air pressure Pa, the manifold temperature Tb, the manifold humidity Hrb, and the manifold pressure Pb. The detail of the humidity detecting EGR rate calculation unit 52 will be described later.

<The Opening Area Learning Value Calculation Unit 53>

The opening area learning value calculation unit 53 is provided with a humidity detecting recirculation flow rate calculation unit 56 that calculates a humidity detecting recirculation flow rate Qeh, which is a flow rate of the recirculation exhaust gas, based on the humidity detecting EGR rate Regr and the intake air flow rate Qa. In Embodiment 1, the humidity detecting EGR rate Regr calculated by the humidity detecting EGR rate calculation unit 52 is an absolute EGR rate which is a ratio of the recirculation exhaust gas to the sum of the intake air and the recirculation exhaust gas, as shown in the equation (11) described later. Then, as shown in an equation (3), the humidity detecting recirculation flow rate calculation unit 56 converts the absolute EGR rate Regr into a relative EGR rate Regrr which is a ratio of the recirculation exhaust gas to the intake air, and then calculates a humidity detecting recirculation flow rate Qeh by multiplying the intake air flow rate Qa to the relative EGR rate Regrr.

$\begin{matrix} {{{{Re}\mspace{11mu} {grr}} = \frac{{Re}\mspace{11mu} {gr}}{1 - {{Re}\mspace{11mu} {gr}}}}{{Qeh} = {{Re}\mspace{11mu} {{grr} \cdot {Qa}}}}} & (3) \end{matrix}$

The opening area learning value calculation unit 53 is provided with a humidity detecting opening area calculation unit 57 that calculates a humidity detecting opening area Segrh which is an opening area of the EGR valve 22 which realizes the humidity detecting recirculation flow rate Qeh, and a learning value calculation unit 58 that calculates a learning value ΔSegrL of the opening area of the EGR valve 22 based on the humidity detecting opening area Segrh.

In Embodiment 1, the humidity detecting opening area calculation unit 57 calculates a sonic velocity Ae of the exhaust gas at the exhaust path 17 side of the EGR valve 22 based on the temperature Tex of the exhaust gas. The humidity detecting opening area calculation unit 57 calculates a density ρe of the exhaust gas at the exhaust path 17 side of the EGR valve 22 based on the temperature Tex of the exhaust gas and the pressure Pex of the exhaust gas. Then, the humidity detecting opening area calculation unit 57 calculates the humidity detecting opening area Segrh based on the manifold pressure Pb, the pressure Pex of the exhaust gas, the sonic velocity Ae of the exhaust gas, the density ρe of the exhaust gas, and the humidity detecting recirculation flow rate Qeh.

Specifically, the humidity detecting opening area calculation unit 57 calculates the humidity detecting opening area Segrh which realizes the humidity detecting recirculation flow rate Qeh, by use of an orifice flow rate calculation equation which is a fluid-mechanics theoretical formula for a compressible fluid, in which the flow in the vicinity of the EGR valve 22 is regarded as flows before and after a throttle valve. The theoretical formula for the flow rate Qe [g/s] of the recirculation exhaust gas that flows through the EGR valve 22, regarded as a throttle valve, is derived as represented in the equation (4), from the energy conservation law, the isentropic flow relational equation, the sonic velocity relational equation, and the state equation.

$\begin{matrix} {{{Qe} = {{{Ae} \cdot {\rho e} \cdot {Segr} \cdot \sigma}\; e}}{{{Ae} = \sqrt{\kappa \cdot R \cdot {Tex}}},{{\rho \; e} = \frac{Pex}{R \cdot {Tex}}}}{{\sigma \; e} = \sqrt{\frac{2}{\kappa - 1}\left\lbrack {\left( \frac{Pb}{Pex} \right)^{\frac{2}{\kappa}} - \left( \frac{Pb}{Pex} \right)^{\frac{\kappa + 1}{\kappa}}} \right\rbrack}}} & (4) \end{matrix}$

Where κ is a specific heat ratio of the recirculation exhaust gas, and a preliminarily set value is used. R is a gas constant of the recirculation exhaust gas, and a preliminarily set value is used. Segr is an opening area of the EGR valve 22. σe is a dimensionless flow rate constant that varies in accordance with a pressure ratio Pb/Pex of the flow at downstream side of (after) the EGR valve 22 to the flow at the upstream side of (before) the EGR valve 22.

The humidity detecting opening area calculation unit 57 calculates a sonic velocity Ae of the exhaust gas based on the temperature Tex of the exhaust gas by use of the second equation of the equation (4). The humidity detecting opening area calculation unit 57 calculates a density ρe of the exhaust gas based on the temperature Tex of the exhaust gas and the pressure Pex of the exhaust gas by use of the third equation of the equation (4).

The humidity detecting opening area calculation unit 57 calculates the dimensionless flow rate constant σe corresponding to the present pressure ratio Pb/Pex of the pressure Pex of the exhaust gas and the manifold pressure Pb, by use of a flow rate constant characteristic data in which the relationship between the pressure ratio Pb/Pex of the pressure Pex of the exhaust gas and the manifold pressure Pb, and the dimensionless flow rate constant σe is preliminarily set based on the fourth equation of the equation (4).

Then, as shown in the equation (5) which is obtained by rearranging the first equation of the equation (4) with regard to the opening area Segr, the humidity detecting opening area calculation unit 57 calculates the humidity detecting opening area Segrh by dividing the humidity detecting recirculation flow rate Qeh by the sonic velocity Ae, the density ρe, and the dimensionless flow rate constant σe of the exhaust gas.

$\begin{matrix} {{Segrh} = \frac{Qeh}{{{Ae} \cdot \rho}\; {e \cdot \sigma}\; e}} & (5) \end{matrix}$

The learning value calculation unit 58 calculates a base opening area Segrb corresponding to the present opening degree Oe of the EGR valve 22, by use of a base opening characteristic data in which the relationship between the base opening area Segrb of the EGR valve 22 and the opening degree Oe of the EGR valve 22 is preliminarily set. Then, the learning value calculation unit 58 calculates a learning value ΔSegrL of opening area, based on the comparison result between the base opening area Segrb and the humidity detecting opening area Segrh.

In Embodiment 1, as shown in the equation (6), the learning value calculation unit 58 calculates a difference ΔSegrh of opening area between the humidity detecting opening area Segrh and the base opening area Segrb; calculates a value obtained by applying an averaging processing (in this example, a first-order-lag filtering processing) to the difference ΔSegrh of opening area as the learning value ΔSegrL of opening area; and stores the learning value ΔSegrL to the storage apparatus 91 such as nonvolatile RAM.

ΔSegrh(n)=Segrh(n)−Segrb(n)

ΔSegrL(n)=Ks·ΔSegrL(n−1)+(1−Ks)·ΔSegrh(n)  (6)

Where (n) denotes the value in the present calculation cycle, and (n−1) denotes the value in the immediately previous calculation cycle. Ks denotes a filter gain in the first-order lag filtering processing and is preliminarily set to a value corresponding to the time constant. Averaging processing, such as a moving-averaging processing, may be performed instead of the first order lag filtering processing, for example. A ratio of opening area and the like may be used instead of the difference ΔSegrh of opening area.

The averaging processing can reduce the influence of the response delay of the manifold humidity sensor 10, the response deviation between the humidity detecting recirculation flow rate Qeh and the flow rate of the recirculation exhaust gas which passes the EGR valve 22 because the humidity detecting EGR rate Regr is an EGR rate after mixing in the intake manifold 12, other disturbance factors, and the like; and stability and accuracy of the learning value ΔSegrL of opening area can be improved.

Alternatively, the learning value calculation unit 58 may increase or decrease the learning value ΔSegrL of opening area, based on the comparison result between the humidity detecting opening area Segrh and the learned opening area SegrL described later. For example, the learning value calculation unit 58 increases the learning value ΔSegrL of opening area, in the case where the humidity detecting opening area Segrh is larger than the learned opening area SegrL; and the learning value calculation unit 58 decreases the learning value ΔSegrL of opening area, in the case where the humidity detecting opening area Segrh is smaller than the learned opening area SegrL.

The learning value calculation unit 58 may calculate the learning value ΔSegrL of opening area for each operating point of the opening degree Oe of the EGR valve 22. For example, the learning value calculation unit 58 stores the learning value ΔSegrL of opening area to the storage apparatus 91, such as nonvolatile RAM, for each opening degree section where the opening degree Oe of the EGR valve 22 was preliminarily divided into a plurality of sections; then the learning value calculation unit 58 reads out the learning value ΔSegrL of the opening degree section corresponding to the present opening degree Oe of the EGR valve 22 from the storage apparatus 91, and updates the learning value ΔSegrL by the difference ΔSegrh of opening area. Accordingly, the same number of the learning value ΔSegrL of opening area as the number of the opening degree section is provided.

The learning value calculation unit 58 permits a update of the learning value ΔSegrL of opening area using the equation (6) in the case of determining that a change of the EGR rate is small and in a steady state; and the learning value calculation unit 58 prohibits the update of the learning value ΔSegrL of opening area using the equation (6) and holds the learning value ΔSegrL of opening area in the case of determining that the change of the EGR rate is large and in a transient state. For example, in the case where a period, in which a change amount of the opening degree Oe of the EGR valve 22 is less than or equal to a preliminarily set EGR determination value and a change amount of the opening degree of the throttle valve 6 is less than or equal to a preliminarily set throttle determination value, passes a preliminarily set decision period, the learning value calculation unit 58 determines that the change of the EGR rate is in the steady state; otherwise, the learning value calculation unit 58 determines that the change of the EGR rate is in the transient state.

Such learning permission conditions can reduce the influence of the response delay of the manifold humidity sensor 10, the response deviation between the humidity detecting recirculation flow rate Qeh and the flow rate of the recirculation exhaust gas which passes the EGR valve 22 because the humidity detecting EGR rate Regr is an EGR rate after mixing in the intake manifold 12, and the like; and accuracy of the learning value ΔSegrL of opening area can be improved.

<The Recirculation Exhaust Gas Calculation Unit for Control 54>

The recirculation exhaust gas calculation unit for control 54 is provided with a learned opening area calculation unit 59 that calculates a learned opening area SegrL of the EGR valve 22 corresponding to the present opening degree Oe of the EGR valve 22 using the learning value ΔSegrL of opening area, and a recirculation flow rate calculation unit 60 for control that calculates a flow rate Qes of the recirculation exhaust gas for control used for control of the internal combustion engine 1 based on the learned opening area SegrL.

Since the humidity detecting recirculation flow rate Qeh which is delayed in a response due to the response delay of the humidity sensor is not used directly, but the learning value ΔSegrL of opening area calculated based on the humidity detecting recirculation flow rate Qeh is used, while suppressing that a response delay due to the humidity sensor causes in the flow rate Qes of the recirculation exhaust gas for control, the calculation accuracy of the flow rate Qes of the recirculation exhaust gas for control can be enhanced.

In Embodiment 1, the learned opening area calculation unit 59 calculates the learned opening area SegrL by correcting the base opening area Segrb by the learning value ΔSegrL of opening area. Then, the recirculation flow rate calculation unit 60 for control calculates the flow rate Qes of the recirculation exhaust gas for control, based on the learned opening area SegrL, the manifold pressure Pb, the pressure Pex of the exhaust gas, the sonic velocity Ae of the exhaust gas, and the density ρe of the exhaust gas.

The learned opening area calculation unit 59 calculates a value obtained by adding the learning value ΔSegrL of opening area to the base opening area Segrb corresponding to the present opening degree Oe of the EGR valve 22, as the learned opening area SegrL, as shown in the equation (7). As is the case with the learning value calculation unit 58, the learned opening area calculation unit 59 calculates the base opening area Segrb corresponding to the present opening degree Oe of the EGR valve 22 by use of the base opening characteristic data described above. The base opening area Segrb which the learning value calculation unit 58 calculated may be used.

SegL=Segrb+ΔSegrL  (7)

In the case where the opening area learning value calculation unit 53 is configured to calculate the learning value ΔSegrL of opening area for each operating point of the opening degree Oe of the EGR valve 22, the learned opening area calculation unit 59 calculates the flow rate Qes of the recirculation exhaust gas for control using the learning value ΔSegrL of opening area corresponding to the present opening degree Oe of the EGR valve 22. Specifically, the learned opening area calculation unit 59 reads out the learning value ΔSegrL of opening degree section corresponding to the present opening degree Oe of the EGR valve 22 from the storage apparatus 91, and calculates a value obtained by adding the read learning value ΔSegrL to the base opening area Segrb, as the learned opening area SegrL.

The recirculation flow rate calculation unit for control 60 calculates the flow rate Qes of the recirculation exhaust gas for control realized by the learned opening area SegrL by use of the orifice flow rate calculation equation of the equation (4), as is the case with the opening area learning value calculation unit 53. As is the case with the humidity detecting opening area calculation unit 57, the recirculation flow rate calculation unit for control 60 calculates the sonic velocity Ae of the exhaust gas based on the temperature Tex of the exhaust gas by use of the second equation of the equation (4). As is the case with the humidity detecting opening area calculation unit 57, the recirculation flow rate calculation unit for control 60 calculates the density ρe of the exhaust gas based on the temperature Tex of the exhaust gas and the pressure Pex of the exhaust gas by use of the third equation of the equation (4). As is the case with the humidity detecting opening area calculation unit 57, the recirculation flow rate calculation unit for control 60 calculates the dimensionless flow rate constant σe corresponding to the present pressure ratio Pb/Pex of the pressure Pex of the exhaust gas and the manifold pressure Pb, by use of the flow rate constant characteristic data described above. The sonic velocity Ae of the exhaust gas, the density ρe of the exhaust gas, and the dimensionless flow rate constant ce, which the humidity detecting opening area calculation unit 57 calculated, may be used.

Then, the recirculation flow rate calculation unit for control 60 calculates a value obtained by multiplying the sonic velocity Ae, the density ρe, and the dimensionless flow rate constant σe of the exhaust gas to the learned opening area SegrL, as the flow rate Qes of the recirculation exhaust gas for control, by use of the equation (8) corresponding to the first equation of the equation (4).

Qes=Ae·ρe·SegrL·σe  (8)

The recirculation exhaust gas calculation unit for control 54 is provided with an EGR rate calculation unit for control 61 that calculates an EGR rate Regrs for control based on the flow rate Qes of the recirculation exhaust gas for control. The EGR rate calculation unit for control 61 calculates a recirculation exhaust gas amount QES [g/stroke] which recirculates to the intake manifold 12 in one stroke period (for example, the interval of BTDC5degCA) based on the flow rate Qes [g/s] of the recirculation exhaust gas for control, as shown in the equation (9); and applies first-order-lag filtering processing, which simulates a delay in the intake manifold 12 (surge tank), to the recirculation exhaust gas amount QES, so as to calculate a cylinder intake recirculation exhaust gas amount QESc [g/stroke] which is an amount of the recirculation exhaust gas taken into the cylinder 25 in one stroke period. The EGR rate calculation unit for control 61 calculates the recirculation exhaust gas amount QES by multiplying one stroke period AT to the flow rate Qes [g/s] of the recirculation exhaust gas for control, for example.

QESc(n)=KCCA·QESc(n−1)+(1−KCCA)·QES(n)

QES(n)=Qes(n)·ΔT(n)  (9)

Where KCCA is the preliminarily set filter gain, and the same value as one of the equation (1) can be used.

The EGR rate calculation unit for control 61 calculates the EGR rate Regrs for control based on the cylinder intake air amount QAc and the cylinder intake recirculation exhaust gas amount QESc. In Embodiment 1, as shown in the equation (10), the EGR rate calculation unit for control 61 calculates, as the EGR rate Regrs for control, a relative EGR rate which divided the cylinder intake recirculation exhaust gas amount QESc by the cylinder intake air amount QAc. An absolute EGR rate may be calculated.

$\begin{matrix} {{{Re}\mspace{11mu} {grs}} = \frac{QESc}{QAc}} & (10) \end{matrix}$

<Recirculation Amount Utilization Control Unit 55>

The recirculation amount utilization control unit 55 controls the internal combustion engine 1 using the recirculation exhaust gas flow rate Qes for control which is calculated by the recirculation exhaust gas calculation unit for control 54. In Embodiment 1, the recirculation amount utilization control unit 55 performs at least one or more of a change of the ignition timing, a change of the opening degree Oe of the EGR valve 22, and a calculation of an output torque of the internal combustion engine 1, based on the flow rate Qes of the recirculation exhaust gas for control.

For example, the recirculation amount utilization control units 55 calculates the ignition timing, based on the rotation speed Ne of the internal combustion engine 1, the charging efficiency Ec, and the EGR rate Regrs for control. The recirculation amount utilization control unit 55 calculates a target EGR rate based on driving condition such as the rotational speed Ne of the internal combustion engine 1 and the charging efficiency Ec, and increases or decreases the opening degree Oe of the EGR valve 22 so that the EGR rate Regrs for control approaches the target EGR rate. By improving the calculation accuracy of the recirculation exhaust gas flow rate Qes for control, the setting precision of the ignition timing and the control accuracy of the EGR rate can be improved, and the control accuracy of the combustion condition, the output torque, the thermal efficiency, the NOx generation amount, and the like of the internal combustion engine 1 can be improved.

The recirculation amount utilization control units 55 calculates the thermal efficiency, based on the rotation speed Ne of the internal combustion engine 1, the charging efficiency Ec, and the EGR rate Regrs for control. Then, the recirculation amount utilization control unit 55 calculates an indicated mean effective pressure by multiplying the thermal efficiency to the calorific value of the fuel supplied to the cylinder 25, and calculates the output torque of the internal combustion engine 1 based on the indicated mean effective pressure. The recirculation amount utilization control unit 55 changes the ignition timing, the intake air amount, and the recirculation exhaust gas amount based on the output torque of the internal combustion engine 1, or transmits the output torque of the internal combustion engine 1 to other controllers such as the transmission controller 81 and makes other controllers use the output torque for a torque control of the whole vehicle.

1-2-1. Detailed Explanation of the Humidity Detecting EGR Rate Calculation Unit 52

Next, the humidity detecting EGR rate calculation unit 52 is explained in detail.

1-2-1-1. Theoretical Derivation of the Calculation Method for the Humidity Detecting EGR Rate Regr

At first, theoretical derivation of a calculation method for the humidity detecting EGR rate Regr will be explained. In Embodiment 1, the humidity detecting EGR rate Regr is the absolute EGR rate, and is the ratio of the exhaust gas (recirculation exhaust gas) recirculated into the intake manifold 12 to the gas in the intake manifold 12. In the following description, the humidity detecting EGR rate Regr also be referred to simply as the EGR rate Regr. By use of the equation (11), the EGR rate Regr can be calculated based on the CO2 concentration.

$\begin{matrix} {{{Re}\mspace{11mu} {gr}} = \frac{{CO}_{2{\_ in}} - {CO}_{2{\_ a}}}{{CO}_{2{\_ ex}} - {CO}_{2{\_ a}}}} & (11) \end{matrix}$

Where CO2_in is the concentration [vol %] of CO2 in the gas inside the intake manifold 12; CO2_ex is the concentration [vol %] of CO2 in the exhaust gas inside the exhaust path 17; CO2_a is the concentration [vol %] of CO2 in the intake air. In general, the concentration of CO2 in the intake air (the atmospheric air) is approximately 0.038 [vol %].

Hereinafter, by paying attention to the respective numbers of moles of molecules and the partial pressures of gases in a combustion chemical reaction formula, a relational equation among the respective numbers of moles of molecules, the partial pressures, the CO2 concentration, and the EGR rate Regr will be derived. The combustion chemical reaction formula for a hydrocarbon at a time when the fuel of the internal combustion engine 1 is gasoline, for example, is expressed by the equation (12).

$\begin{matrix} {\left. {{C_{n}H_{m}} + {\left( {n + \frac{m}{4}} \right) \cdot O_{2}}}\rightarrow{{n \cdot {CO}_{2}} + {{\frac{m}{2} \cdot H_{2}}O}} \right.~} & (12) \end{matrix}$

In the case where when it is assumed that the average molecular formula for gasoline is C7H14 and that the composition of air is “oxygen (O2): nitrogen (N2)=21:79”, the gasoline and the air combust together at the theoretical air-fuel ratio AF0, the combustion chemical reaction formula is expressed by the equation (13). In the equation (13), the respective numbers of moles of carbon dioxide (CO2) and water vapor (H2O), which are produced by the combustion, are the same, i.e., 14.

2.C₇H₁₄+21.O₂+79.N₂→14.CO₂+14.H₂O+79.N₂  (13)

However, actual intake air includes carbon dioxide (CO2) and water vapor (H2O); thus, when it is assumed that the respective numbers of moles of carbon dioxide (CO2) and water vapor (H2O) are α and β, the combustion chemical reaction formula is expressed by the equation (14). The equation (14) is a pure combustion chemical reaction formula in which no recirculation exhaust gas is taken into consideration.

2.C₇H₁₄+21.O₂+79.N₂+α.CO₂+β.H₂O→α(14+α).CO₂+(14+β).H₂O+79.N₂  (14)

In the following analysis of the number of moles, the number of moles of gasoline in the left-hand side of the equation (14) is approximated with “0” because it is small in comparison with the number of total moles. Accordingly, the pre-combustion gas represented in the left-hand side of the equation (14) becomes equal to the intake air. Because the number of total moles of the intake air represented in the left-hand side of the equation (14) is (100+α+β) and the number of total moles of the exhaust gas represented in the right-hand side of the equation (14) is (107+α+β), the respective numbers of moles at the left-hand side and at the right-hand side are different from each other, strictly speaking; however, in this example, it is assumed that both the respective numbers of total moles at the left-hand side and at the right-hand side are (M+α+β).

FIG. 4 represents the state of partial pressures of gasses in the intake manifold 12 at a time when exhaust gas is recirculated to the intake manifold 12. The gas in the intake manifold 12 is a mixture gas including intake air taken from the atmosphere and recirculation exhaust gas recirculated through the EGR path 21; the partial pressure of the intake air in the manifold pressure is expressed by “P_new”, and the partial pressure of the recirculation exhaust gas is expressed by “P_egr”.

As represented in the left-hand side of the equation (14), the intake air is composed of nitrogen (N2), oxygen (O2), carbon dioxide (CO2), and water vapor (H2O). Strictly speaking, although included in the intake air, other substances are neglected because their content is extremely small. Pvn denotes the partial pressure of water vapor included in the intake air.

As represented in the right-hand side of the equation (14), the recirculation exhaust gas is composed of nitrogen (N2), carbon dioxide (CO2), and water vapor (H2O). Strictly speaking, although included in the intake air, other substances are neglected because their content is extremely small. The carbon dioxide (CO2) includes carbon dioxide (CO2) that is produced through combustion and carbon dioxide (CO2) that has originally been included in the intake air; the water vapor (H2O) includes water vapor (H2O) that is produced through combustion and water vapor (H2O) that has originally been included in the intake air. Accordingly, the partial pressure of the water vapor produced through combustion is expressed by Pve, and the partial pressure of the water vapor that has originally been included in the intake air is expressed by Pvr.

When the respective CO2 concentrations in the equation (11) for calculating the EGR rate Regr are expressed by the respective ratios of the partial pressure P_new of the intake air and the partial pressure P_egr of the recirculation exhaust gas to the manifold pressure Pb, represented in FIG. 4, and the mole fraction of CO2 in the intake air represented in the left-hand side of the equation (14) or in the exhaust gas represented in the right-hand side of the equation (14), the equation (15) is yielded. Specifically, as represented in FIG. 4, the concentration CO2_in of CO2 in the gas inside the intake manifold 12 is the total of the concentration of CO2 in the intake air inside the intake manifold 12, the concentration of CO2 that is produced through combustion and included in the recirculation exhaust gas, and the concentration of CO2 included in the intake air. The concentration of CO2 in the intake air inside the intake manifold 12 is obtained by multiplying the ratio of the partial pressure P_new of the intake air to the manifold pressure Pb by the mole fraction (c/(M+α+β)) of CO2 in the intake air, represented in the left-hand side of the equation (14). The concentration of CO2 produced through combustion and CO2 in the intake air that are included in the recirculation exhaust gas is obtained by multiplying the ratio of the partial pressure P_egr of the recirculation exhaust gas to the manifold pressure Pb by the mole fraction ((14+α)/(M+α+β)) of CO2 in the exhaust gas, represented in the right-hand side of the equation (14). The concentration CO2_ex of CO2 in the exhaust gas is the mole fraction ((14+α)/(M+α+β)) of CO2 in the exhaust gas, represented in the right-hand side of the equation (14). The concentration CO2_a of CO2 in the intake air is the mole fraction (α/(M+α+β)) of CO2 in the intake air, represented in the left-hand side of the equation (14).

$\begin{matrix} {{{CO}_{2{\_ in}} = {{\frac{P\_ new}{Pb} \cdot \frac{\alpha}{M + \alpha + \beta}} + {\frac{P\_ egr}{Pb} \cdot \frac{14 + \alpha}{M + \alpha + \beta}}}}{{CO}_{2{\_ ex}} = \frac{14 + \alpha}{M + \alpha + \beta}}{{CO}_{2{\_ a}} = \frac{\alpha}{M + \alpha + \beta}}} & (15) \end{matrix}$

When the respective concentrations of CO2 in the equation (15) are substituted for the equation (11) and the equation (1) is rearranged, the EGR rate Regr is given by the ratio (P_egr/Pb) of the partial pressure P_egr of the recirculation exhaust gas to the manifold pressure Pb, as represented in the equation (16).

$\begin{matrix} \begin{matrix} {{{Re}\; {gr}} = {\frac{{CO}_{2{\_ in}} - {CO}_{2{\_ a}}}{{CO}_{2{\_ ex}} - {CO}_{2{\_ a}}} = \frac{\begin{matrix} {\begin{pmatrix} {{\frac{P\_ new}{Pb} \cdot \frac{\alpha}{M + \alpha + \beta}} +} \\ {\frac{P\_ egr}{Pb} \cdot \frac{14 + \alpha}{M + \alpha + \beta}} \end{pmatrix} -} \\ \left( \frac{\alpha}{M + \alpha + \beta} \right) \end{matrix}}{\left( \frac{14 + \alpha}{M + \alpha + \beta} \right) - \left( \frac{\alpha}{M + \alpha + \beta} \right)}}} \\ {= \frac{P\_ egr}{Pb}} \end{matrix} & (16) \end{matrix}$

Next, a relational equation will be derived by paying attention to the mole fractions of the gases. By use of the mole fraction (β/(M+α+β)) of water vapor in the intake air represented in the left-hand side of the equation (14), the mole fraction χva of the water vapor in the intake air can be expressed as in the equation (17).

$\begin{matrix} {{\chi \; {va}} = \frac{\beta}{M + \alpha + \beta}} & (17) \end{matrix}$

As represented in FIG. 4 and the equation (18), the mole fraction χvb of the water vapor in the gas inside the intake manifold 12 is the total of the mole fraction (the first item in the right-hand side of the equation (18)) of the water vapor included in the intake air inside the intake manifold 12 and the mole fraction (the second item in the right-hand side of the equation (18)) of the water vapor produced through combustion and the water vapor in the intake air, which are included in the recirculation exhaust gas. As represented in the first item of the right-hand side of the equation (18), the mole fraction of the water vapor included in the intake air inside the intake manifold 12 is obtained by multiplying the ratio (P_new/Pb) of the partial pressure P_new of the intake air to the manifold pressure Pb by the mole fraction (β/(M+α+β)) of the water vapor in the intake air, represented in the left-hand side of the equation (14). As represented in the second item of the right-hand side of the equation (18), the mole fraction of the water vapor included in the intake air and the water vapor produced through combustion, which are included in the recirculation exhaust gas, is obtained by multiplying the ratio (P_egr/Pb) of the partial pressure P_egr of the recirculation exhaust gas to the manifold pressure Pb by the mole fraction ((14+β/(M+α+β)) of the water vapor in the exhaust gas, represented in the right-hand side of the equation (14). Accordingly, the mole fraction χvb of the water vapor in the intake manifold 12 can be expressed as in the equation (18), by use of these mole fractions of water vapor.

$\begin{matrix} \begin{matrix} {{\chi \; {vb}} = {{\frac{P\_ new}{Pb} \cdot \frac{\beta}{M + \alpha + \beta}} + {\frac{P\_ egr}{Pb} \cdot \frac{14 + \alpha}{M + \alpha + \beta}}}} \\ {= {\frac{\beta}{M + \alpha + \beta} + {\frac{P\_ egr}{Pb} \cdot \frac{14}{M + \alpha + \beta}}}} \end{matrix} & (18) \end{matrix}$

The rearranged first item (β/(M+α+β)) of the right-hand side of the equation (18) is equal to the mole fraction χva of the water vapor in the intake air, represented in the equation (17); thus, when the equation (17) is substituted for the equation (18) and then the equation (18) is rearranged, the equation (19) is obtained. As represented in the equation (19), the subtraction mole fraction (χvb−χva) obtained by subtracting the mole fraction χva of the water vapor in the intake air from the mole fraction χvb of the water vapor in the intake manifold 12 is equal to the value obtained by multiplying the ratio (P_egr/Pb) of the partial pressure P_egr of the recirculation exhaust gas to the manifold pressure Pb by the mole fraction χvex (=14/(M+α+β)) of inner-exhaust-gas combustion-produced water vapor, which is the mole fraction of water vapor produced through combustion in the exhaust gas, obtained from the right-hand side of the equation (14). The subtraction mole fraction (χvb−χva) becomes equal to the mole fraction χvegr of the water vapor, produced through combustion, in the recirculation exhaust gas (referred to as the mole fraction χvegr of combustion-produced water vapor in the intake manifold 12). Thus, the equation (19) suggests that the mole fraction χvegr of combustion-produced water vapor in the intake manifold 12 is equal to a multiplication value obtained by multiplying the ratio (P_egr/Pb) of the partial pressure P_egr of the recirculation exhaust gas to the manifold pressure Pb by the mole fraction χvex of inner-exhaust-gas combustion-produced water vapor.

${{\chi \; {vb}} - {\chi \; {va}}} = {{\frac{P\_ egr}{Pb} \cdot \chi}\; {vex}\begin{matrix} {{{{\,\chi}\; {vex}} = \frac{14}{M + \alpha + \beta}}{{\chi \; {vegr}} = {{\chi \; {vb}} - {\chi \; {va}}}}} & (19) \end{matrix}}$

From the derivation result of the equation (15), (P_egr/Pb) in the equation (19) is equal to the EGR rate Regr; thus, when the equation (15) is substituted for the equation (19) and then the equation (19) is rearranged, the equation (20) is obtained. Accordingly, the EGR rate Regr becomes a value obtained by dividing the mole fraction χvegr of combustion-produced water vapor in the intake manifold 12, calculated by subtracting the mole fraction χva of the water vapor in the intake air from the mole fraction χvb of the water vapor in the intake manifold 12, by the mole fraction χvex of inner-exhaust-gas combustion-produced water vapor. In other words, it is suggested that by dividing the mole fraction χvegr of combustion-produced water vapor in the intake manifold 12 by the mole fraction χvex of inner-exhaust-gas combustion-produced water vapor, the total mole fraction of all the recirculation exhaust gas in the gas inside the intake manifold 12 is obtained and that the mole fraction of the recirculation exhaust gas is equal to the ratio (P_egr/Pb) of the partial pressure P_egr of the recirculation exhaust gas to the manifold pressure Pb, i.e., the recirculation exhaust gas.

$\begin{matrix} {{{{Re}\; {gr}} = {\left( {{\chi \; {vb}} - {\chi \; {va}}} \right) \cdot \frac{1}{\chi \; {vex}}}}{{\chi \; {vex}} = \frac{14}{M + \alpha + \beta}}{{\chi \; {vegr}} = {{\chi \; {vb}} - {\chi \; {va}}}}} & (20) \end{matrix}$

By rearranging the equation (17) with respect to β, the number β of moles of water vapor in the intake air can be expressed by the equation (21).

$\begin{matrix} {\beta = {\frac{\chi \; {va}}{1 - {\chi \; {va}}} \cdot \left( {M + \alpha} \right)}} & (21) \end{matrix}$

By substituting the equation (21) for the equation (20) and then rearranging the equation (20), the equation (22) is obtained. In the equation (22), the number α of moles of CO2 in the intake air is set to 0.038, which is a nominal value.

$\begin{matrix} {\mspace{79mu} {{{{Re}\; {gr}} = {\left( {{\chi \; {vb}} - {\chi \; {va}}} \right) \cdot \frac{1}{\chi \; {vex}}}}{{\chi \; {vex}} = {{\frac{14}{M + \alpha + {\frac{\chi \; {va}}{1 - {\chi \; {va}}} \cdot \left( {M + \alpha} \right)}} \cdot \left( {1 - {\chi \; {va}}} \right)} = {{\frac{14}{M + \alpha} \cdot \left( {1 - {\chi \; {va}}} \right)} = {\frac{14}{107 + 0.038} \cdot \left( {1 - {\chi \; {va}}} \right)}}}}}} & (22) \end{matrix}$

Thus, it is conceivable that from the derivation result of the equation (22), the EGR rate Regr can be calculated based on the mole fraction χvb of water vapor in the intake manifold 12 and the mole fraction χva of water vapor in intake air.

As represented in the equation (23), the mole fraction χvb of water vapor in the intake manifold 12 becomes theoretically equal to the ratio (Pvb/Pb) of the partial pressure Pvb of inner-manifold water vapor, which is the partial pressure of water vapor included in the gas inside the intake manifold 12, to the manifold pressure Pb. The mole fraction χva of water vapor in the intake air becomes theoretically equal to the ratio (Pva/Pa) of the partial pressure Pva of inner-intake-air water vapor, which is the partial pressure of water vapor included in the intake air, to the intake-air pressure Pa.

$\begin{matrix} {{{\chi \; {vb}} = \frac{Pvb}{Pb}}{{\chi \; {va}} = \frac{Pva}{Pa}}} & (23) \end{matrix}$

Thus, it is conceivable that as represented in the equation (24) obtained by substituting the equation (23) for the equation (22), the EGR rate Regr can be calculated by detecting the inner-manifold water vapor partial pressure ratio (Pvb/Pb), which is the ratio of the partial pressure Pvb of the inner-manifold water vapor to the manifold pressure Pb, and the inner-intake-air water vapor partial pressure ratio (Pva/Pa), which is the ratio of the partial pressure Pva of the inner-intake-air water vapor to the intake-air pressure Pa.

$\begin{matrix} {{{Re}\; {gr}} = {\left( {\frac{Pvb}{Pb} - \frac{Pva}{Pa}} \right) \cdot \frac{M + \alpha}{14} \cdot \frac{1}{1 - \frac{Pva}{Pa}}}} & (24) \end{matrix}$

1-2-1-2. Configuration of the Humidity Detecting EGR Rate Calculation Unit 52

The humidity detecting EGR rate calculation unit 52 calculates the humidity detecting EGR rate Regr based on the intake-air temperature Ta, the intake-air humidity Hra, the intake-air pressure Pa, the manifold temperature Tb, the manifold humidity Hrb, and the manifold pressure Pb. In Embodiment 1, the humidity detecting EGR rate calculation unit 52 is provided with an inner-manifold water vapor ratio calculation unit 70, an inner-intake-air water vapor ratio calculation unit 71, and a final EGR rate calculation unit 72, as shown in FIG. 2.

Based on the manifold humidity Hrb and the manifold temperature Tb which are detected by the driving-condition detection unit 51, the inner-manifold water vapor ratio calculation unit 70 calculates the partial pressure Pvb of the inner-manifold water vapor which is the partial pressure of water vapor included in the gas inside the intake manifold 12, so as to calculate the inner-manifold water vapor partial pressure ratio (Pvb/Pb) which is the ratio of the partial pressure Pvb of the inner-manifold water vapor to the manifold pressure Pb. Based on the intake-air humidity Hra and the intake-air temperature Ta which are detected by the driving-condition detection unit 51, the inner-intake-air water vapor ratio calculation unit 71 calculates the partial pressure Pva of inner-intake-air water vapor, which is the partial pressure of water vapor included in the intake air, so as to calculate the inner-intake-air water vapor partial pressure ratio (Pva/Pa), which is the ratio of the partial pressure Pva of the inner-intake-air water vapor to the intake-air pressure Pa.

Then, based on the inner-manifold water vapor partial pressure ratio (Pvb/Pb) and the inner-intake-air water vapor partial pressure ratio (Pva/Pa), the final EGR rate calculation unit 72 calculates the humidity detecting EGR rate Regr, which is the ratio of recirculation exhaust gas, recirculated into the intake manifold 12, to the gas in the intake manifold 12.

This configuration makes it possible that based on the detection values of the pressure, the temperature, and the humidity of each of gas and intake air in the intake manifold 12, the inner-manifold water vapor partial pressure ratio (Pvb/Pb) and the inner-intake-air water vapor partial pressure ratio (Pva/Pa) are calculated and that based on these two water vapor partial pressure ratios, the humidity detecting EGR rate Regr is accurately calculated. The details will be explained below.

<The Inner-Manifold Water Vapor Ratio Calculation Unit 70>

As described above, based on the manifold humidity Hrb and the manifold temperature Tb, the inner-manifold water vapor ratio calculation unit 70 calculates the partial pressure Pvb of the inner-manifold water vapor, which is the partial pressure of water vapor included in the gas inside the intake manifold 12, so as to calculate the inner-manifold water vapor partial pressure ratio (Pvb/Pb), which is the ratio of the partial pressure Pvb of the inner-manifold water vapor to the manifold pressure Pb.

As described above, the inner-manifold water vapor partial pressure ratio (Pvb/Pb) is theoretically equal to the water vapor mole fraction χvb, which is the ratio of the number of moles of water vapor to the number of total moles of gas in the intake manifold 12. Therefore, the inner-manifold water vapor partial pressure ratio (Pvb/Pb) can be referred to also as the water vapor mole fraction χvb.

In Embodiment 1, the inner-manifold water vapor ratio calculation unit 70 calculates a saturated water vapor pressure Psb in the intake manifold 12, based on the manifold temperature Tb, and then calculates the partial pressure Pvb of the inner-manifold water vapor by multiplying the saturated water vapor pressure Psb by the manifold humidity Hrb, which is a relative humidity.

Specifically, the inner-manifold water vapor ratio calculation unit 70 calculates the saturated water vapor pressure Psb in the intake manifold 12 by use of the Tetens equation represented in the equation (25).

$\begin{matrix} {{Psb} = {6.1078 \times 10^{(\frac{7.5 \times {Tb}}{{Tb} + 237.3})}}} & (25) \end{matrix}$

The inner-manifold water vapor ratio calculation unit 70 may calculate the saturated water vapor pressure Psb corresponding to the manifold temperature Tb, by use of a characteristic data in which the relationship between the temperature and the saturated water vapor pressure is preliminarily set based on the equation (25).

Then, as represented in the equation (26), the inner-manifold water vapor ratio calculation unit 70 calculates the partial pressure Pvb of the inner-manifold water vapor, by multiplying the saturated water vapor pressure Psb by the manifold humidity Hrb [% RH], and then calculates the inner-manifold water vapor partial pressure ratio (Pvb/Pb), by dividing the partial pressure Pvb of the inner-manifold water vapor by the manifold pressure Pb.

$\begin{matrix} {{Pvb} = {{Psb} \cdot \frac{Hrb}{100}}} & (26) \end{matrix}$

<The Inner-Intake-Air Water Vapor Ratio Calculation Unit 71>

As described above, based on the intake-air humidity Hra and the intake-air temperature Ta, the inner-intake-air water vapor ratio calculation unit 71 calculates the partial pressure Pva of inner-intake-air water vapor, which is the partial pressure of water vapor included in the intake air, so as to calculate the inner-intake-air water vapor partial pressure ratio (Pva/Pa), which is the ratio of the partial pressure Pva of the inner-intake-air water vapor to the intake-air pressure Pa.

As described above, the inner-intake-air water vapor partial pressure ratio (Pva/Pa) is theoretically equal to the water vapor mole fraction χva, which is the ratio of the number of moles of water vapor to the number of total moles of intake air. Thus, the inner-intake-air water vapor partial pressure ratio (Pva/Pa) can be referred to also as the inner-intake-air water vapor mole fraction χva.

In Embodiment 1, based on the intake-air temperature Ta, the inner-intake-air water vapor ratio calculation unit 71 calculates a saturated water vapor pressure Psa of intake air and then calculates the inner-intake-air water vapor partial pressure Pva by multiplying the saturated water vapor pressure Psa by the intake-air humidity Hra, which is a relative humidity.

Specifically, the inner-intake-air water vapor ratio calculation unit 71 calculates the saturated water vapor pressure Psa of the intake air by use of the Tetens equation represented in the equation (27).

$\begin{matrix} {{Psa} = {6.1078 \times 10^{(\frac{7.5 \times {Ta}}{{Ta} + 237.3})}}} & (27) \end{matrix}$

The inner-intake-air water vapor ratio calculation unit 71 may calculate the saturated water vapor pressure Psa corresponding to the intake-air temperature Ta by use of a characteristic data in which the relationship between the temperature and the saturated water vapor pressure is preliminarily set.

Then, as represented in the equation (28), the inner-intake-air water vapor ratio calculation unit 71 calculates the inner-intake-air water vapor partial pressure Pva by multiplying the saturated water vapor pressure Psa by the intake-air humidity Hra [% RH], and then calculates the inner-intake-air water vapor partial pressure ratio (Pva/Pa), by dividing the partial pressure Pva of the inner-intake-air water vapor by the intake-air pressure Pa.

$\begin{matrix} {{Pva} = {{Psa} \cdot \frac{Hra}{100}}} & (28) \end{matrix}$

<The Final EGR Rate Calculation Unit 72>

As described above, based on the inner-manifold water vapor partial pressure ratio (Pvb/Pb) and the inner-intake-air water vapor partial pressure ratio (Pva/Pa), the final EGR rate calculation unit 72 calculates the humidity detecting EGR rate Regr, which is the ratio of recirculation exhaust gas, recirculated into the intake manifold 12, to the gas in the intake manifold 12.

The final EGR rate calculation unit 72 calculates the humidity detecting EGR rate Regr by use of the equation (29) based on the equation (24) above. In other words, the final EGR rate calculation unit 72 calculates a subtraction partial pressure ratio by subtracting the inner-intake-air water vapor partial pressure ratio (Pva/Pa) from the inner-manifold water vapor partial pressure ratio (Pvb/Pb), multiplies the subtraction partial pressure ratio by a preliminarily set conversion constant Kr, then calculates, as the humidity detecting EGR rate Regr, a value by dividing the multiplication value by a subtraction value obtained by subtracting the inner-intake-air water vapor partial pressure ratio (Pva/Pa) from “1”.

$\begin{matrix} {{{Regr} = {\left( {\frac{Pvb}{Pb} - \frac{Pva}{Pa}} \right) \cdot {Kr} \cdot \frac{1}{1 - \frac{Pva}{Pa}}}}{{Kr} = {\frac{M + \alpha}{14} = \frac{107 + 0.038}{14}}}} & (29) \end{matrix}$

Based on the respective numbers of moles of molecules in the chemical reaction formula at a time when a fuel and wet air combust together, the conversion constant Kr is preliminarily set as represented in the equation (29). Specifically, the conversion constant Kr is preliminarily set to a fixed value obtained by dividing a value (107+0.038, in this example), calculated by subtracting the number β of moles of inner-intake-air water vapor from the number (M+α+β) of total moles of exhaust gas in the right-hand side (combusted gas) of the combustion chemical reaction formula represented in the equation (14), by the number of moles of the combustion-produced water vapor (14, in this example). The conversion constant Kr may be a value other than the value represented in the equation (29), for example, a value obtained through adjustment based on an experimental value. Because sufficiently smaller than M, the number α of moles of inner-intake-air carbon dioxide may be regarded as zero.

The calculation by the equation (29) will be expressed based on physical quantities. The equation (30) is obtained by modifying the equation (29). As represented in the equation (30), the final EGR rate calculation unit 72 calculates an inner-manifold combustion-produced water vapor partial pressure ratio (Pvegr/Pb), which is the ratio of the partial pressure Pvegr (referred to as an inner-intake-manifold combustion-produced water vapor partial pressure Pvegr) of combustion-produced water vapor included in recirculation exhaust gas to the manifold pressure Pb, by subtracting the inner-intake-air water vapor partial pressure ratio (Pva/Pa) from the inner-manifold water vapor partial pressure ratio (Pvb/Pb). Based on the number of moles of molecules in the chemical reaction formula at a time when a fuel and wet air combust together and the inner-intake-air water vapor partial pressure ratio (Pva/Pa), the final EGR rate calculation unit 72 calculates the inner-exhaust-gas combustion-produced water vapor mole fraction χvex, which is the mole fraction of combustion-produced water vapor in exhaust gas. Then, the final EGR rate calculation unit 72 calculates, as the humidity detecting EGR rate Regr, a value by dividing the inner-manifold combustion-produced water vapor partial pressure ratio (Pvegr/Pb) by the inner-exhaust-gas combustion-produced water vapor mole fraction χvex.

$\begin{matrix} {{{Regr} = {\frac{Pvegr}{Pb} \cdot \frac{1}{\chi \; {vex}}}}{\frac{Pvegr}{Pb} = {\frac{Pvb}{Pb} - \frac{Pva}{Pa}}}{{\chi \; {vex}} = {{Kr}\; {2 \cdot \left( {1 - \frac{Pva}{Pa}} \right)}}}{{{Kr}\; 2} = {\frac{14}{M + \alpha} = \frac{14}{107 + 0.038}}}} & (30) \end{matrix}$

Speaking in detail, the final EGR rate calculation unit 72 calculates, as the inner-exhaust-gas combustion-produced water vapor mole fraction χvex, a value by multiplying a preliminarily set mole conversion constant Kr2 by a subtraction value calculated by subtracting the inner-intake-air water vapor partial pressure ratio (Pva/Pa) from “1”. As is the case with the foregoing conversion constant Kr, the mole conversion constant Kr2 is preliminarily set based on an experimental value, the respective numbers of moles of molecules in the chemical reaction formula at a time when a fuel and wet air combust together, or the like. The calculation equation for the inner-exhaust-gas combustion-produced water vapor mole fraction χvex is derived by substituting the equation (21), derived based on the combustion chemical reaction formula represented in the equation (14), for the equation (20), derived based on the combustion chemical reaction formula represented in the equation (14).

As described above, by dividing the inner-manifold combustion-produced water vapor partial pressure ratio (Pvegr/Pb) by the inner-exhaust-gas combustion-produced water vapor mole fraction χvex, the ratio (P_egr/Pb) of the partial pressure P_egr of recirculation exhaust gas to the manifold pressure Pb is obtained and hence the humidity detecting EGR rate Regr is obtained.

1-2-2. Flowchart

The procedure (the control method for the internal combustion engine 1) of the processing by the controller 50 according to Embodiment 1 will be explained based on the flowchart represented in FIG. 5. The processing represented in the flowchart in FIG. 5 is recurrently implemented, for example, every constant operation cycle while the computing processing unit 90 implements software (a program) stored in the storage apparatus 91.

In the step S01, the driving-condition detection unit 51 implements driving condition detection processing (a driving condition detection step) for, as mentioned above, detecting the driving condition of the internal combustion engine 1. The driving-condition detection unit 51 detects the manifold pressure Pb, the manifold temperature Tb, the manifold humidity Hrb, the intake-air pressure Pa, the intake-air temperature Ta, the intake-air humidity Hra, the intake air flow rate Qa, the opening degree Oe of the EGR valve 22, and the like.

Next, in the step S02, the humidity detecting EGR rate calculation unit 52 implements a humidity detecting EGR rate calculation processing (a humidity detecting EGR rate calculation Step) for, as mentioned above, calculating the humidity detecting EGR rate Regr based on the intake-air temperature Ta, the intake-air humidity Hra, the intake-air pressure Pa, the manifold temperature Tb, the manifold humidity Hrb, and the manifold pressure Pb. In the step S02, an inner-manifold water vapor ratio calculation processing (an inner-manifold water vapor ratio calculation step) performed by the inner-manifold water vapor ratio calculation unit 70, an inner-intake-air water vapor ratio calculation processing (an inner-intake-air water vapor ratio calculation step) performed by the inner-intake-air water vapor ratio calculation unit 71, and a final EGR rate calculation processing (a final EGR rate calculation Step) performed by the final EGR rate calculation unit 72 are performed in order.

In the step S03, the opening area learning value calculation unit 53 implements an opening area learning value calculation processing (an opening area learning value calculation step) for, as mentioned above, calculating the humidity detecting recirculation flow rate Qeh based on the humidity detecting EGR rate Regr and the intake air flow rate Qa, calculating the humidity detecting opening area Segrh, which is an opening area of the EGR valve 22 which realizes the humidity detecting recirculation flow rate Qeh, and calculating the learning value ΔSegrL of the opening area of the EGR valve 22 based on the humidity detecting opening area Segrh.

In the step S04, the recirculation exhaust gas calculation unit for control 54 implements a recirculation exhaust gas calculation processing for control (a recirculation exhaust gas calculation step for control) for, as mentioned above, calculating the learned opening area SegrL of the EGR valve 22 corresponding to the present opening degree Oe of the EGR valve 22 using the learning value ΔSegrL of opening area, and calculating the flow rate Qes of the recirculation exhaust gas for control used for control of the internal combustion engine 1 based on the learned opening area SegrL.

In the step S05, the recirculation amount utilization control unit 55 implements a recirculation amount utilization control processing (a recirculation amount utilization control step) for, as mentioned above, controlling the internal combustion engine 1 using the recirculation exhaust gas flow rate Qes for control.

2. Embodiment 2

Next, the controller 50 according to Embodiment 2 will be explained. The explanation for constituent parts the same as those in Embodiment 1 will be omitted. The respective basic configurations and processing methods of the internal combustion engine 1 and the controller 50 according to Embodiment 2 are the same as those of the internal combustion engine 1 and the controller 50 according to Embodiment 1; however, Embodiment 2 is different from Embodiment 1 in that the humidity detecting EGR rate calculation unit 52 changes the humidity detecting EGR rate Regr in accordance with the air-fuel ratio AF of the internal combustion engine 1.

2-1. Extension of the Calculation Method for the Humidity Detecting EGR Rate Regr to the Case of Rich- or Lean-Air-Fuel Ratio

The calculation method for the humidity detecting EGR rate Regr according to Embodiment 1 is based on the combustion chemical reaction formula represented in the equation (14) at a time when a fuel and wet air combust together at the theoretical air-fuel ratio AF0. Hereinafter, the equation derivation will be extended to the case where the air-fuel ratio of wet air to a fuel is leaner than the theoretical air-fuel ratio AF0 or to the case where the air-fuel ratio is richer that the theoretical air-fuel ratio AF0.

An excess air ratio λ is the ratio of the air-fuel ratio AF to the theoretical air-fuel ratio AF0, as represented in the equation (31). When λ=1, the air-fuel ratio AF is equal to the theoretical air-fuel ratio AF0; when λ<1, the air-fuel ratio AF is rich; when λ>1, the air-fuel ratio AF is lean.

$\begin{matrix} {\lambda = \frac{AF}{{AF}\; 0}} & (31) \end{matrix}$

<In the Case of Rich Air-Fuel Ratio>

The equation (32) represents the combustion chemical reaction formula at a time when the air-fuel ratio AF is rich (λ<1).

2.C₇H₁₄+λ·{21.O₂+79.N₂+α.CO₂+β.H₂O}→λ·{(14+α).CO₂+(14+β).H₂O+79.N₂}+(1−λ).2.C₇H₁₄  (32)

In this situation, as represented by the last term in the right-hand side of the equation (32), it is assumed that in the case of rich combustion, uncombusted gasoline is directly exhausted with its original molecules. In practice, it is conceivable that due to the combustion temperature in the cylinder 25, uncombusted gasoline is decomposed into methane (CH4), ethane (C2H6), and the like that each have a molecular weight smaller than that of gasoline; however, because the volume concentration of gasoline is small and the decomposition of gasoline may not provide a substantial effect, the decomposition of gasoline is not taken into consideration in this embodiment.

As represented in the equation (33), the concentration CO2_ex of CO2 in exhaust gas is equal to the ratio (the mole fraction of CO2) of the number of moles of CO2 to the number of total moles in the exhaust gas represented in the right-hand side of the equation (32). Because being small in comparison with the number of total moles, the number ((1−λ)·2) of moles of uncombusted gasoline is approximated with zero. Similarly, the concentration CO2_in of CO2 in the gas inside the intake manifold 12 and the concentration CO2_a of CO2 in the intake air are also obtained. From the equation (33), the concentration of CO2 at a time when the air-fuel ratio is rich becomes equal to the concentration of CO2 represented in the equation (15) at a time of the theoretical air-fuel ratio AF0. Therefore, in the case where the air-fuel ratio is rich, the EGR rate Regr becomes the ratio (P_egr/Pb) of the partial pressure P_egr of the recirculation exhaust gas to the manifold pressure Pb, as is the case with the equation (16).

$\begin{matrix} {\begin{matrix} {{CO}_{2{\_ ex}} = \frac{\lambda \cdot \left( {14 + \alpha} \right)}{{\lambda \cdot \left\{ {M + \alpha + \beta} \right\}} + {\left( {1 - \lambda} \right) \cdot 2}}} \\ {\cong \frac{\lambda \cdot \left( {14 + \alpha} \right)}{\lambda \cdot \left\{ {M + \alpha + \beta} \right\}}} \\ {= \frac{14 + \alpha}{M + \alpha + \beta}} \end{matrix}\begin{matrix} {{CO}_{2{\_ {in}}} = {{\frac{P\_ new}{Pb} \cdot \frac{\lambda \cdot \alpha}{\lambda \cdot \left\{ {M + \alpha + \beta} \right\}}} + {\frac{P\_ egr}{Pb} \cdot}}} \\ {\frac{\lambda \cdot \left( {14 + \alpha} \right)}{{\lambda \cdot \left\{ {M + \alpha + \beta} \right\}} + {\left( {1 - \lambda} \right) \cdot 2}}} \\ {\cong {{\frac{P\_ new}{Pb} \cdot \frac{\alpha}{M + \alpha + \beta}} + {\frac{P\_ egr}{Pb} \cdot \frac{14 + \alpha}{M + \alpha + \beta}}}} \end{matrix}{{\begin{matrix} {{CO}_{2{\_ a}} = \frac{\lambda \cdot \alpha}{\lambda \cdot \left\{ {M + \alpha + \beta} \right\}}} \\ {= \frac{\alpha}{M + \alpha + \beta}} \end{matrix}\therefore{Regr}} = \frac{P\_ egr}{Pb}}} & (33) \end{matrix}$

In the case where the air-fuel ratio is rich, the mole fraction χva of water vapor in the intake air and the mole fraction χvb of water vapor in the gas in the intake manifold 12 are represented by the equation (34) by use of the mole fraction of water vapor in the left-hand side or the right-hand side of the equation (32) and the like, as is the case with the respective mole fractions at a time of the theoretical air-fuel ratio AF0 represented in the equations (17) and (18).

$\begin{matrix} {\begin{matrix} {{\chi \; {va}} = \frac{\lambda \cdot \beta}{\lambda \cdot \left\{ {M + \alpha + \beta} \right\}}} \\ {= \frac{\beta}{M + \alpha + \beta}} \end{matrix}\begin{matrix} {{\chi \; {vb}} = {{\frac{P\_ new}{Pb} \cdot \frac{\lambda \cdot \beta}{\lambda \cdot \left\{ {M + \alpha + \beta} \right\}}} + {\frac{P\_ egr}{Pb} \cdot}}} \\ {\frac{\lambda \cdot \left( {14 + \beta} \right)}{{\lambda \cdot \left\{ {M + \alpha + \beta} \right\}} + {\left( {1 - \lambda} \right) \cdot 2}}} \\ {\cong {{\frac{P\_ new}{Pb} \cdot \frac{\beta}{M + \alpha + \beta}} + {\frac{P\_ egr}{Pb} \cdot \frac{14 + \beta}{M + \alpha + \beta}}}} \\ {= {\frac{\beta}{M + \alpha + \beta} + {\frac{P\_ egr}{Pb} \cdot \frac{14}{M + \alpha + \beta}}}} \end{matrix}} & (34) \end{matrix}$

From the equation (34), the mole fractions χva and χvb of water vapor at a time when the air-fuel ratio is rich become equal to those represented in the equations (17) and (18), respectively, at a time of the theoretical air-fuel ratio AF0. When at a time of a rich air-fuel ratio (λ<1), uncombusted gasoline is produced and hence the numerator, i.e., the number of moles of combustion-produced water vapor is obtained by multiplying the number of moles of combustion-produced water vapor at a time of the theoretical air-fuel ratio AF0, i.e., 14 by λ; however, because the respective numbers of moles of the molecules in the intake air are multiplied by λ, the denominator, i.e., the number of total moles of exhaust gas is also obtained by multiplying the number of total moles of exhaust gas at a time of the theoretical air-fuel ratio AF0 by λ. As a result, λ in the numerator and λ in the denominator cancel each other and hence the respective values the same as those at a time of the theoretical air-fuel ratio AF0 can be obtained. Accordingly, when the air fuel ratio is rich, the EGR rate Regr can be calculated by the equation (34), which is the equation for the case where the air fuel ratio is the theoretical air-fuel ratio AF0.

<In the Case of Lean Air-Fuel Ratio>

The equation (35) represents the combustion chemical reaction formula at a time when the air-fuel ratio AF is lean (λ>1).

2.C₇H₁₄+λ.{21.O₂+79.N₂+α.CO₂+β.H₂O}→(14+λ·α).CO₂+(14+λ·β).H₂O+λ·79.N₂+21·(λ−1).O₂  (35)

As represented in the equation (36), the concentration CO2_ex of CO2 in exhaust gas is equal to the ratio (the mole fraction of CO2) of the number of moles of CO2 to the number of total moles in the exhaust gas represented in the right-hand side of the equation (35). Because being small in comparison with the number of total moles, the number of moles “7” is approximated with zero. Similarly, the concentration CO2 in of CO2 in the gas inside the intake manifold 12 and the concentration CO2_a of CO2 in the intake air are also obtained. From the equation (36), the concentration of CO2 at a time when the air-fuel ratio is lean differs from the concentration of CO2 represented in the equation (15) at a time of the theoretical air-fuel ratio AF0.

$\begin{matrix} {\begin{matrix} {{CO}_{2{\_ ex}} = \frac{\left( {14 + {\lambda \cdot \alpha}} \right)}{\left( {14 + {\lambda \cdot \alpha}} \right) + \left( {14 + {\lambda \cdot \beta}} \right) + {\lambda \cdot 79} + {21 \cdot \left( {\lambda - 1} \right)}}} \\ {= \frac{14 + {\lambda \cdot \alpha}}{{\lambda \cdot \left( {M + \alpha + \beta} \right)} + 7}} \\ {\cong \frac{14 + {\lambda \cdot \alpha}}{\lambda \cdot \left( {M + \alpha + \beta} \right)}} \end{matrix}\begin{matrix} {{CO}_{2{\_ {in}}} = {{\frac{P\_ new}{Pb} \cdot \frac{\lambda \cdot \alpha}{\lambda \cdot \left\{ {M + \alpha + \beta} \right\}}} + {\frac{P\_ egr}{Pb} \cdot}}} \\ {\frac{14 + {\lambda \cdot \alpha}}{{\lambda \cdot \left( {M + \alpha + \beta} \right)} + 7}} \\ {\cong {{\frac{P\_ new}{Pb} \cdot \frac{\alpha}{M + \alpha + \beta}} + {\frac{P\_ egr}{Pb} \cdot \frac{14 + {\lambda \cdot \alpha}}{\lambda \cdot \left( {M + \alpha + \beta} \right)}}}} \end{matrix}\begin{matrix} {{CO}_{2{\_ a}} = \frac{\lambda \cdot \alpha}{\lambda \cdot \left\{ {M + \alpha + \beta} \right\}}} \\ {= \frac{\alpha}{M + \alpha + \beta}} \end{matrix}} & (36) \end{matrix}$

Then, by substituting the respective concentrations of CO2 in the equation (36) for the equation (11) and then rearranging the equation (1), the equation (37) is obtained. In the case where the air-fuel ratio is lean, the EGR rate Regr also becomes the ratio (P_egr/Pb) of the partial pressure P_egr of the recirculation exhaust gas to the manifold pressure Pb, as is the case with the equation (16).

$\begin{matrix} \begin{matrix} {{Regr} = \frac{{CO}_{2{\_ in}} - {CO}_{2{\_ a}}}{{CO}_{2{\_ ex}} - {CO}_{2{\_ a}}}} \\ {= \frac{\begin{matrix} {\left( {{\frac{P\_ new}{Pb} \cdot \frac{\alpha}{M + \alpha + \beta}} + {\frac{P\_ egr}{Pb} \cdot \frac{14 + {\lambda \cdot \alpha}}{\lambda \cdot \left( {M + \alpha + \beta} \right)}}} \right) -} \\ \left( \frac{\alpha}{M + \alpha + \beta} \right) \end{matrix}}{\left( \frac{14 + {\lambda\_\alpha}}{\lambda \cdot \left( {M + \alpha + \beta} \right)} \right) - \left( \frac{\alpha}{M + \alpha + \beta} \right)}} \\ {= \frac{P\_ egr}{Pb}} \end{matrix} & (37) \end{matrix}$

In the case where the air-fuel ratio is lean, the mole fraction χva of water vapor in the intake air and the mole fraction χvb of water vapor in the gas in the intake manifold 12 are represented by the equation (38) by use of the mole fraction of water vapor in the left-hand side or the right-hand side of the equation (35) and the like, as is the case with the respective mole fractions at a time of the theoretical air-fuel ratio AF0 represented in the equations (17) and (18).

$\begin{matrix} {\begin{matrix} {{\chi \; {va}} = \frac{\lambda \cdot \beta}{\lambda \cdot \left\{ {M + \alpha + \beta} \right\}}} \\ {= \frac{\beta}{M + \alpha + \beta}} \end{matrix}\begin{matrix} {{\chi \; {vb}} = {{\frac{P\_ new}{Pb} \cdot \frac{\lambda \cdot \beta}{\lambda \cdot \left\{ {M + \alpha + \beta} \right\}}} + {\frac{P\_ egr}{Pb} \cdot}}} \\ {\frac{14 + {\lambda \cdot \beta}}{{\lambda \cdot \left\{ {M + \alpha + \beta} \right\}} + 7}} \\ {\cong {{\frac{P\_ new}{Pb} \cdot \frac{\beta}{M + \alpha + \beta}} + {\frac{P\_ egr}{Pb} \cdot \frac{14 + {\lambda \cdot \beta}}{\lambda \cdot \left( {M + \alpha + \beta} \right)}}}} \\ {= {\frac{\beta}{M + \alpha + \beta} + {\frac{P\_ egr}{Pb} \cdot \frac{14}{\lambda \cdot \left( {M + \alpha + \beta} \right)}}}} \end{matrix}} & (38) \end{matrix}$

By rearranging the equation (38), as is the case with the equation (19), the equation (39) is obtained. The inner-exhaust-gas combustion-produced water vapor mole fraction χvexL at a time when the air-fuel ratio is lean becomes a value obtained by dividing the inner-exhaust-gas combustion-produced water vapor mole fraction χvex, represented in the equation (19), at a time of the theoretical air-fuel ratio AF0 by the excess air ratio λ. The reason for that is the following: when at a time of a lean air-fuel ratio (λ>1), the fuel completely combusts and hence the numerator, i.e., the number of moles of combustion-produced water vapor is the same as the number of moles of combustion-produced water vapor at a time of the theoretical air-fuel ratio AF0, i.e., 14; in contrast, because the respective numbers of moles of the molecules in the intake air are multiplied by λ, the denominator, i.e., the number of total moles of exhaust gas is obtained by multiplying the number of total moles of exhaust gas at a time of the theoretical air-fuel ratio AF0 by λ.

$\begin{matrix} {{{{\chi \; {vb}} - {\chi \; {va}}} = {{\frac{P\_ egr}{Pb} \cdot \chi}\; {vexL}}}{{\chi \; {vexL}} = {\frac{1}{\lambda} \cdot \frac{14}{M + \alpha + \beta}}}{{\chi \; {vegr}} = {{\chi \; {vb}} - {\chi \; {va}}}}} & (39) \end{matrix}$

By substituting the equation (37) for the equation (39) and then rearranging the equation (39), the equation (40) is obtained.

$\begin{matrix} {{{Regr} = {\left( {{\chi \; {vb}} - {\chi \; {va}}} \right) \cdot \frac{1}{\chi \; {vexL}}}}{{\chi \; {vexL}} = {{{\frac{1}{\lambda} \cdot \frac{14}{M + \alpha + \beta}}\chi \; {vegr}} = {{\chi \; {vb}} - {\chi \; {va}}}}}} & (40) \end{matrix}$

The number β of moles of water vapor in the intake air is given by the equation (21), as is the case with the theoretical air-fuel ratio AF0; thus, by substituting the equation (21) for the equation (40) and rearranging the equation (40), the equation (41) is obtained.

$\begin{matrix} {{{Regr} = {\left( {{\chi \; {vb}} - {\chi \; {va}}} \right) \cdot \frac{1}{\chi \; {vexL}}}}\begin{matrix} {{\chi \; {vexL}} = {\frac{1}{\lambda} \cdot \frac{14}{M + \alpha + {\frac{\chi \; {va}}{1 - {\chi \; {va}}} \cdot \left( {M + \alpha} \right)}} \cdot \left( {1 - {\chi \; {va}}} \right)}} \\ {= {\frac{1}{\lambda} \cdot \frac{14}{M + \alpha} \cdot \left( {1 - {\chi \; {va}}} \right)}} \end{matrix}} & (41) \end{matrix}$

Then, by substituting the equation (23) for the equation (41) and then rearranging the equation (41), the equation (42) is obtained. Therefore, in the case of a lean air-fuel ratio, the EGR rate Regr can be calculated by multiplying the EGR rate Regr, represented in the equation (24), at a time of the theoretical air-fuel ratio AF0 by the excess air ratio λ.

$\begin{matrix} {{Regr} = {\lambda \cdot \left( {\frac{Pvb}{Pb} - \frac{Pva}{Pa}} \right) \cdot \frac{M + \alpha}{14} \cdot \frac{1}{1 - \frac{Pva}{Pa}}}} & (42) \end{matrix}$

The foregoing derivation results will be summarized in the equation (43). In the case where the air-fuel ratio AF of the internal combustion engine 1 is the theoretical air-fuel ratio (AF=AF0) or rich (AF<AF0), the EGR rate Regr can be calculated, through the equation (24), based on the inner-manifold water vapor partial pressure ratio (Pvb/Pb) and the inner-intake-air water vapor partial pressure ratio (Pva/Pa). In the case where the air-fuel ratio AF of the internal combustion engine 1 is lean (AF>AF0), the EGR rate Regr can be calculated by, as represented in the equation (42), further multiplying the EGR rate Regr calculated through the equation (24) by the excess air ratio λ.

$\begin{matrix} {\left. 1 \right)\mspace{14mu} {in}\mspace{14mu} {the}\mspace{14mu} {case}\mspace{14mu} {of}\mspace{14mu} {theoretical}\mspace{14mu} {air}\text{/}{fuel}\mspace{14mu} {ratio}\mspace{14mu} {or}\mspace{14mu} {rich}} & \; \\ {{\left. {{{Regr} = {\left( {\frac{Pvb}{Pb} - \frac{Pva}{Pa}} \right) \cdot \frac{M + \alpha}{14} \cdot \frac{1}{1 - \frac{Pva}{Pa}}}}2} \right)\mspace{14mu} {in}\mspace{14mu} {the}\mspace{14mu} {case}\mspace{14mu} {of}\mspace{14mu} {lean}}{{Regr} = {\lambda \cdot \left( {\frac{Pvb}{Pb} - \frac{Pva}{Pa}} \right) \cdot \frac{M + \alpha}{14} \cdot \frac{1}{1 - \frac{Pva}{Pa}}}}} & (43) \end{matrix}$

2-2. The Configuration of the Final EGR Rate Calculation Unit 72

Thus, the final EGR rate calculation unit 72 according to Embodiment 2 changes the humidity detecting EGR rate Regr in accordance with the air-fuel ratio AF of the internal combustion engine 1.

As is the case with the equation (29) or (30) in Embodiment 1, the final EGR rate calculation unit 72 calculates the humidity detecting EGR rate Regr0 at a time when it is assumed that the air-fuel ratio AF of the internal combustion engine 1 is the theoretical air-fuel ratio, based on the inner-manifold water vapor partial pressure ratio (Pvb/Pb) and the inner-intake-air water vapor partial pressure ratio (Pva/Pa). As represented in the equation (44), in the case where the air-fuel ratio AF is leaner than the theoretical air-fuel ratio AF0, the final EGR rate calculation unit 72 calculates, as the final humidity detecting EGR rate Regr, a value by multiplying the humidity detecting EGR rate Regr0 at a time when it is assumed that the air-fuel ratio is the theoretical air-fuel ratio AF0 by the excess air ratio λ obtained by dividing the air-fuel ratio AF by the theoretical air-fuel ratio AF0. On the other hand, in the case where the air-fuel ratio AF is richer than the theoretical air-fuel ratio AF0, the final EGR rate calculation unit 72 directly calculates, as the final humidity detecting EGR rate Regr, the humidity detecting EGR rate Regr0 at a time when it is assumed that the air-fuel ratio is the theoretical air-fuel ratio AF0.

1) in the case of lean

Regr=λ·Regr0

λ=AF/AF0  (44)

2) in the case of theoretical air/fuel ratio or rich

Regr=Regr0

Expression will be made based on physical quantities. As represented in the equation (45), the final EGR rate calculation unit 72 calculates, as is the case with the equation (40) in Embodiment 1, the inner-exhaust-gas combustion-produced water vapor mole fraction χvex0 at a time when it is assumed that the air-fuel ratio AF of the internal combustion engine 1 is the theoretical air-fuel ratio, based on the respective numbers of moles of molecules in the chemical reaction formula at a time when the fuel and wet air combust together at the theoretical air-fuel ratio AF0 and the inner-intake-air water vapor partial pressure ratio (Pva/Pa). The final EGR rate calculation unit 72 calculates the inner-manifold combustion-produced water vapor partial pressure ratio (Pvegr/Pb) by subtracting the inner-intake-air water vapor partial pressure ratio (Pva/Pa) from the inner-manifold water vapor partial pressure ratio (Pvb/Pb).

$\begin{matrix} {{{{\chi \; {vex}\; 0} = {{Kr}\; {2 \cdot \left( {1 - \frac{Pva}{Pa}} \right)}}},{{{Kr}\; 2} = {\frac{14}{M + \alpha} = \frac{14}{107 + 0.038}}}}{\frac{Pvegr}{Pb} = {\frac{Pvb}{Pb} - \frac{Pva}{Pa}}}} & (45) \end{matrix}$

As represented in the equation (46), in the case where the air-fuel ratio AF is leaner than the theoretical air-fuel ratio AF0, the final EGR rate calculation unit 72 calculates, as the inner-exhaust-gas combustion-produced water vapor mole fraction χvexL at a time of a lean air-fuel ratio, a value by dividing the inner-exhaust-gas combustion-produced water vapor mole fraction χvex0 at a time when it is assumed that the air-fuel ratio is the theoretical air-fuel ratio AF0 by the excess air ratio A. Then, the final EGR rate calculation unit 72 calculates, as the humidity detecting EGR rate Regr, a value by dividing the inner-manifold combustion-produced water vapor partial pressure ratio (Pvegr/Pb) by the inner-exhaust-gas combustion-produced water vapor mole fraction χvexL at a time of a lean air-fuel ratio. On the other hand, in the case where the air-fuel ratio AF is richer than the theoretical air-fuel ratio AF0, the final EGR rate calculation unit 72 directly adopts, as the inner-exhaust-gas combustion-produced water vapor mole fraction χvex0 at a time of a rich air-fuel ratio, the inner-exhaust-gas combustion-produced water vapor mole fraction χvex0 at a time when it is assumed that the air-fuel ratio is the theoretical air-fuel ratio AF0. Then, the final EGR rate calculation unit 72 calculates, as the humidity detecting EGR rate Regr, a value by dividing the inner-manifold combustion-produced water vapor partial pressure ratio (Pvegr/Pb) by the inner-exhaust-gas combustion-produced water vapor mole fraction χvex0 at a time of a rich air-fuel ratio.

$\begin{matrix} {\left. 1 \right)\mspace{14mu} {in}\mspace{14mu} {the}\mspace{14mu} {case}\mspace{14mu} {of}\mspace{14mu} {lean}} & \; \\ {{\left. {{{\chi \; {vexL}} = {{\frac{1}{\lambda} \cdot \chi}\; {vex}\; 0}}{{Regr} = {\frac{Pvegr}{Pb} \cdot \frac{1}{\chi \; {vexL}}}}2} \right)\mspace{14mu} {in}\mspace{14mu} {the}\mspace{14mu} {case}\mspace{14mu} {of}\mspace{14mu} {theoretical}\mspace{14mu} {air}\text{/}{fuel}\mspace{14mu} {ratio}\mspace{14mu} {or}\mspace{14mu} {rich}}{{Regr} = {\frac{Pvegr}{Pb} \cdot \frac{1}{\chi \; {vex}\; 0}}}} & (46) \end{matrix}$

Based on the air-fuel ratio AF detected by the air-fuel ratio sensor 18, the final EGR rate calculation unit 72 determines whether the air-fuel ratio AF is rich, lean, or the theoretical air-fuel ratio, and calculates the excess air ratio A. Specifically, in the case where the detection value of the air-fuel ratio AF is smaller than a preliminarily set value (e.g., 14.7) of the theoretical air-fuel ratio AF0 (AF<AF0), the final EGR rate calculation unit 72 determines that the air-fuel ratio is rich; in the case where the detection value of the air-fuel ratio AF is larger than the setting value of the theoretical air-fuel ratio AF0 (AF>AF0), the final EGR rate calculation unit 72 determines that the air-fuel ratio is lean; in the case where the detection value of the air-fuel ratio AF is equal to the setting value of the theoretical air-fuel ratio AF0 (AF=AF0), the final EGR rate calculation unit 72 determines that the air-fuel ratio is the theoretical air-fuel ratio. In addition, the final EGR rate calculation unit 72 calculates the excess air ratio λ (=AF/AF0) by dividing the detection value of the air-fuel ratio AF by the preliminarily set value of the theoretical air-fuel ratio AF0.

Alternatively, it may be allowed that based on a fuel correction coefficient Kaf to be utilized in calculating the fuel injection amount, the final EGR rate calculation unit 72 determines whether the air-fuel ratio AF is rich, lean, or the theoretical air-fuel ratio, and calculates the excess air ratio λ. For example, the fuel correction coefficient Kaf is a correction coefficient by which a basic fuel injection amount calculated so as to realize the theoretical air-fuel ratio AF0 is multiplied. In the case where Kaf=1, the final EGR rate calculation unit 72 determines that the air-fuel ratio is the theoretical air-fuel ratio AF0; in the case where Kaf>1, the final EGR rate calculation unit 72 determines that the air-fuel ratio is rich; in the case where Kaf<1, the final EGR rate calculation unit 72 determines that the air-fuel ratio is lean. As the excess air ratio λ (=1/Kaf), the final EGR rate calculation unit 72 calculates the reciprocal of the fuel correction coefficient Kaf.

Alternatively, it may be allowed that based on a control mode for the air-fuel ratio, the final EGR rate calculation unit 72 determines whether the air-fuel ratio AF is rich, lean, or the theoretical air-fuel ratio, and sets the excess air ratio λ. The control modes for the air-fuel ratio include a theoretical air-fuel ratio control mode, a rich control mode, and a lean control mode. In the theoretical air-fuel ratio control mode, in order to raise the purification performance of the three-way catalyst, the air-fuel ratio AF detected by the air-fuel ratio sensor 18 is feedback-controlled so as to be in the vicinity of the theoretical air-fuel ratio AF0. In the rich control mode, which is set at a time when high-load driving is performed, the air-fuel ratio AF is controlled so as to be rich. In the lean control mode, which is set, for example, in order to raise the gasoline mileage, the air-fuel ratio AF is controlled so as to be lean. In the case where the control mode for the air-fuel ratio is the theoretical air-fuel ratio control mode, the rich control mode, or the lean control mode, the final EGR rate calculation unit 72 determines that the air-fuel ratio is the theoretical air-fuel ratio, rich, or lean, as the case may be. As described above, the final EGR rate calculation unit 72 calculates the excess air ratio A, based on the detection value of the air-fuel ratio AF, detected by the air-fuel ratio sensor 18, or the fuel correction coefficient Kaf.

3. Embodiment 3

Next, the controller 50 according to Embodiment 3 will be explained. The explanation for constituent parts the same as those in each of Embodiments 1 and 2 will be omitted. The respective basic configurations and processing methods of the internal combustion engine 1 and the controller 50 according to Embodiment 3 are the same as those of the internal combustion engine 1 and the controller 50 according to each of Embodiments 1 and 2; however, Embodiment 3 is different from each of Embodiments 1 and 2 in that the intake-air humidity sensor 5 is not provided in the internal combustion engine 1 and is also different in the detection methods for the intake-air pressure Pa, the intake-air temperature Ta, and the intake-air humidity Hra.

In each of Embodiments 1 and 2, there has been explained the case where the driving-condition detection unit 51 detects the intake-air pressure Pa, based on the output signal of the intake-air pressure sensor 2, detects the intake-air temperature Ta, based on the output signal of the intake-air temperature sensor 4, and detects the intake-air humidity Hra, based on the output signal of the intake-air humidity sensor 5.

However, in Embodiment 3, as the intake-air pressure Pa, the intake-air humidity Hra, and the intake-air temperature Ta, the driving-condition detection unit 51 detects the manifold pressure Pb, the manifold humidity Hrb, and the manifold temperature Tb, respectively, under the condition that the EGR valve 22 for opening closing the EGR path 21 is closed and hence no exhaust gas is recirculated to the intake manifold 12.

In Embodiment 3, in the case where the EGR valve 22 has continuously been closed for a preliminarily set determination period, the driving-condition detection unit 51 determines that an intake air detection condition has been satisfied. The determination period is set to a period in which after the EGR valve 22 is closed, recirculation exhaust gas in the intake manifold 12 sufficiently decreases. The determination period may be shortened as the intake air amount increases. While the intake air detection condition is satisfied, the driving-condition detection unit 51 detects the intake-air pressure Pa, based on the output signal of the manifold pressure sensor 8, detects the intake-air temperature Ta, based on the output signal of the manifold temperature sensor 9, and detects the intake-air humidity Hra, based on the output signal of the manifold humidity sensor 10. While the EGR valve 22 is closed, no exhaust gas is recirculated into the intake manifold 12; thus, only intake air exists therein. As a result, by use of the pressure, the humidity, and the temperature of the gas in the intake manifold 12, that are detected under the foregoing condition, the inner-intake-air water vapor partial pressure Pva and the inner-intake-air water vapor partial pressure ratio (Pva/Pa) can be calculated.

As is the case with Embodiment 1, based on the intake-air humidity Hra and the intake-air temperature Ta that are detected at a time when the intake air detection condition is satisfied, the inner-intake-air water vapor ratio calculation unit 71 calculates the partial pressure Pva of inner-intake-air water vapor, which is the partial pressure of water vapor included in the intake air, so as to calculate the inner-intake-air water vapor partial pressure ratio (Pva/Pa), which is the ratio of the partial pressure Pva of the inner-intake-air water vapor to the intake-air pressure Pa. When the intake air detection condition is not satisfied, the inner-intake-air water vapor ratio calculation unit 71 holds and outputs the inner-intake-air water vapor partial pressure ratio (Pva/Pa) that has been calculated at a time when the intake air detection condition was satisfied. The humidity of atmospheric air changes more gently than the humidity in the intake manifold 12 changes due to a change in the EGR rate; therefore, even when the held value is utilized, a large estimation error is prevented from occurring in the EGR rate.

OTHER EMBODIMENTS

Lastly, other embodiments of the present invention will be explained. Each of the configurations of embodiments to be explained below is not limited to be separately utilized but can be utilized in combination with the configurations of other embodiments as long as no discrepancy occurs.

(1) In each of the foregoing embodiments, there has been explained the case in which considering the case where as the fuel for the internal combustion engine 1, gasoline is utilized, it is assumed that the average molecular formula of gasoline and the composition of air are given by each of the equations (14), (32), and (35). However, embodiments of the present invention are not limited to the foregoing case. In other words, it may be allowed that in the average molecular formula of gasoline and the composition of air, stricter values are utilized, that the chemical reaction formula and the respective numbers of moles of molecules in the chemical reaction formula are changed, and that the setting values for calculating the humidity detecting EGR rate Regr such as the conversion constant Kr and the mole conversion constant Kr2 are changed. Moreover, it may be allowed that as the fuel of the internal combustion engine 1, for example, light oil, alcohol, natural gas, or the like is utilized, that the average molecular formula of the fuel is changed in accordance with the kind of the fuel, that the chemical reaction formula and the respective numbers of moles of molecules in the chemical reaction formula are changed, and that the setting values for calculating the humidity detecting EGR rate Regr such as the conversion constant Kr and the mole conversion constant Kr2 are changed.

(2) In each of the foregoing embodiments, there has been explained the case in which the driving-condition detection unit 51 detects the intake-air temperature Ta, based on the output signal of the intake-air temperature sensor 4, and detects the intake-air humidity Hra, based on the output signal of the intake-air humidity sensor 5. However, embodiments of the present invention are not limited to the foregoing case. That is to say, it may be allowed that the driving-condition detection unit 51 obtains information on the intake-air humidity Hra and the intake-air temperature Ta from the air conditioner controller 80. The air conditioner controller 80 is a controller for an air conditioner that performs air-conditioning of a vehicle interior and is connected with the controller 50 for the internal combustion engine 1 through a communication wire. The air conditioner controller 80 is provided with a humidity sensor that detects the humidity of atmospheric air to be taken into the air conditioner and a temperature sensor that detects the temperature of atmospheric air, detects the atmospheric-air humidity and the atmospheric-air temperature based on the output signals of the humidity sensor and the temperature sensor, and transmits information on the atmospheric-air humidity and the atmospheric-air temperature to the controller 50.

(3) In each of the foregoing embodiments, there has been explained the case in which the recirculation amount utilization control unit 55 performs, as mentioned above, at least one or more of a change of the ignition timing, a change of the opening degree Oe of the EGR valve 22, and a calculation of an output torque of the internal combustion engine 1, based on the flow rate Qes of the recirculation exhaust gas for control. However, embodiments of the present invention are not limited to the foregoing case. That is to say, it may be allowed that the recirculation amount utilization control unit 55 uses the recirculation exhaust gas flow rate Qes for control of other than these controls, for example, a control of the intake air amount, a control for changing the valve opening and closing timing of one or both of the intake valve 14 and the exhaust valve 15, and the like.

Various modifications and alterations of this invention will be apparent to those skilled in the art without departing from the scope and spirit of this invention, and it should be understood that this is not limited to the illustrative embodiments set forth herein. 

What is claimed is:
 1. A controller for an internal combustion engine that is provided with an intake path and an exhaust path, a throttle valve for opening and closing the intake path, an EGR path for recirculating exhaust gas from the exhaust path to an intake manifold that is part of the intake path and is at the downstream side of the throttle valve, and an EGR valve for opening and closing the EGR path, the controller for the internal combustion engine comprising: a driving-condition detector that detects a manifold pressure, which is a pressure of gas in the intake manifold, a manifold temperature, which is a temperature of gas in the intake manifold, a manifold humidity, which is a humidity of gas in the intake manifold, an intake-air pressure, which is a pressure of intake air to be taken into the intake path, an intake-air temperature, which is a temperature of the intake air, an intake-air humidity, which is a humidity of the intake air, an intake air flow rate, which is a flow rate of the intake air, and an opening degree of the EGR valve; a humidity detecting EGR rate calculator that calculates a humidity detecting EGR rate, which is a ratio of a recirculation exhaust gas, which is the exhaust gas recirculated into the intake manifold, to the intake air, based on the intake-air temperature, the intake-air humidity, the intake-air pressure, the manifold temperature, the manifold humidity, and the manifold pressure; an opening area learning value calculator that calculates a humidity detecting recirculation flow rate which is a flow rate of the recirculation exhaust gas based on the humidity detecting EGR rate and the intake air flow rate, calculates a humidity detecting opening area, which is an opening area of the EGR valve which realizes the humidity detecting recirculation flow rate, and calculates a learning value of the opening area of the EGR valve based on the humidity detecting opening area; and a recirculation exhaust gas calculator for control that calculates a learned opening area of the EGR valve corresponding to the present opening degree of the EGR valve using the learning value of the opening area, and calculates a flow rate of the recirculation exhaust gas for control used for controlling the internal combustion engine based on the learned opening area.
 2. The controller for the internal combustion engine according to claim 1, wherein the driving-condition detector detects a temperature of the exhaust gas at the exhaust path side of the EGR valve, and a pressure of the exhaust gas at the exhaust path side of the EGR valve, wherein the opening area learning value calculator calculates a sonic velocity of the exhaust gas at the exhaust path side of the EGR valve based on the temperature of the exhaust gas, calculates a density of the exhaust gas at the exhaust path side of the EGR valve based on the temperature of the exhaust gas and the pressure of the exhaust gas, and calculates the humidity detecting opening area based on the manifold pressure, the pressure of the exhaust gas, the sonic velocity of the exhaust gas, the density of the exhaust gas, and the humidity detecting recirculation flow rate.
 3. The controller for the internal combustion engine according to claim 2, wherein the opening area learning value calculator calculates a base opening area corresponding to the present opening degree of the EGR valve, by use of a base opening characteristic data in which the relationship between the base opening area of the EGR valve and the opening degree of the EGR valve is preliminarily set, and calculates the learning value of the opening area based on the comparison result between the base opening area and the humidity detecting opening area, wherein the recirculation exhaust gas calculator for control calculates the learned opening area by correcting the base opening area with the learning value of the opening area, and calculates the flow rate of the recirculation exhaust gas for control based on the learned opening area, the manifold pressure, the pressure of the exhaust gas, the sonic velocity of the exhaust gas, and the density of the exhaust gas.
 4. The controller for the internal combustion engine according to claim 1, wherein the opening area learning value calculator calculates the learning value of the opening area for each operating point of the opening degree of the EGR valve, the recirculation exhaust gas calculator for control calculates the flow rate of the recirculation exhaust gas for control using the learning value of the opening area corresponding to the present opening degree of the EGR valve.
 5. The controller for the internal combustion engine according to claim 1, wherein the opening area learning value calculator permits a update of the learning value of the opening area in the case of determining that a change of the EGR rate is small and in a steady state, and prohibits the update of the learning value of the opening area and holds the learning value of the opening area in the case of determining that the change of the EGR rate is large and in a transient state.
 6. The controller for the internal combustion engine according to claim 1, further comprising: a recirculation amount utilization controller performs at least one or more of a calculation of an output torque of the internal combustion engine, a change of an ignition timing, and a change of the opening degree of the EGR valve, based on the flow rate of the recirculation exhaust gas for control.
 7. The controller for the internal combustion engine according to claim 1, the humidity detecting EGR rate calculator includes: an inner-manifold water vapor ratio calculator that calculates an inner-manifold water vapor partial pressure, which is a partial pressure of water vapor included in gas inside the intake manifold, based on the manifold humidity and the manifold temperature, and calculates an inner-manifold water vapor partial pressure ratio, which is the ratio of the inner-manifold water vapor partial pressure to the manifold pressure; an inner-intake-air water vapor ratio calculator that calculates an inner-intake-air water vapor partial pressure, which is a partial pressure of water vapor included in the intake air, based on the intake-air humidity and the intake-air temperature, and calculates an inner-intake-air water vapor partial pressure ratio, which is the ratio of the inner-intake-air water vapor partial pressure to the intake-air pressure; and a final EGR rate calculator that calculates the humidity detecting EGR rate based on the inner-manifold water vapor partial pressure ratio and the inner-intake-air water vapor partial pressure ratio.
 8. The controller for the internal combustion engine according to claim 7, wherein the final EGR rate calculator calculates a subtraction partial pressure ratio by subtracting the inner-intake-air water vapor partial pressure ratio from the inner-manifold water vapor partial pressure ratio, multiplies the subtraction partial pressure ratio by a preliminarily set conversion constant, and then calculates, as the humidity detecting EGR rate, a value by dividing the multiplication value by a subtraction value obtained by subtracting the inner-intake-air water vapor partial pressure ratio from “1”.
 9. The controller for the internal combustion engine according to claim 7, wherein the final EGR rate calculator calculates an inner-manifold combustion-produced water vapor partial pressure ratio, which is a ratio of the partial pressure of combustion-produced water vapor included in the recirculated exhaust gas to the manifold pressure, by subtracting the inner-intake-air water vapor partial pressure ratio from the inner-manifold water vapor partial pressure ratio; based on the respective numbers of moles of molecules in a chemical reaction formula at a time when a fuel and wet air combust together and the inner-intake-air water vapor partial pressure ratio, the final EGR rate calculator calculates an inner-exhaust-gas combustion-produced water vapor mole fraction, which is a mole fraction of combustion-produced water vapor in the exhaust gas; then, the final EGR rate calculator calculates, as the humidity detecting EGR rate, a value by dividing the inner-manifold combustion-produced water vapor partial pressure ratio by the inner-exhaust-gas combustion-produced water vapor mole fraction.
 10. The controller for the internal combustion engine according to claim 7, wherein the final EGR rate calculator changes the humidity detecting EGR rate in accordance with an air-fuel ratio of the internal combustion engine.
 11. The controller for the internal combustion engine according to claim 7, wherein the final EGR rate calculator calculates the humidity detecting EGR rate at a time when it is assumed that an air-fuel ratio of the internal combustion engine is a theoretical air-fuel ratio, based on the inner-manifold water vapor partial pressure ratio and the inner-intake-air water vapor partial pressure ratio; in the case where the air-fuel ratio is leaner than the theoretical air-fuel ratio, the final EGR rate calculator calculates, as the final humidity detecting EGR rate, a value by multiplying the humidity detecting EGR rate at a time when it is assumed that the air-fuel ratio is the theoretical air-fuel ratio by an excess air ratio obtained by dividing the air-fuel ratio by the theoretical air-fuel ratio; in the case where the air-fuel ratio is richer than the theoretical air-fuel ratio, the final EGR rate calculator directly calculates, as the final humidity detecting EGR rate, the humidity detecting EGR rate at a time when it is assumed that the air-fuel ratio is the theoretical air-fuel ratio.
 12. The controller for the internal combustion engine according to claim 1, wherein as the intake-air pressure, the intake-air humidity, and the intake-air temperature, the driving-condition detector detects the manifold pressure, the manifold humidity, and the manifold temperature under the condition that the EGR valve for the path is closed and hence no exhaust gas is recirculated to the intake manifold.
 13. The controller for the internal combustion engine according to claim 1, wherein the driving-condition detector obtains information on the intake-air humidity and the intake-air temperature from an air conditioner controller.
 14. A control method for an internal combustion engine that is provided with an intake path and an exhaust path, a throttle valve for opening and closing the intake path, an EGR path for recirculating exhaust gas from the exhaust path to an intake manifold that is part of the intake path and is at the downstream side of the throttle valve, and an EGR valve for opening and closing the EGR path, the control method comprising: a driving-condition detecting that detects a manifold pressure, which is a pressure of gas in the intake manifold, a manifold temperature, which is a temperature of gas in the intake manifold, a manifold humidity, which is a humidity of gas in the intake manifold, an intake-air pressure, which is a pressure of intake air to be taken into the intake path, an intake-air temperature, which is a temperature of the intake air, an intake-air humidity, which is a humidity of the intake air, an intake air flow rate, which is a flow rate of the intake air, and an opening degree of the EGR valve; a humidity detecting EGR rate calculating that calculates a humidity detecting EGR rate, which is a ratio of a recirculation exhaust gas, which is the exhaust gas recirculated into the intake manifold, to the intake air, based on the intake-air temperature, the intake-air humidity, the intake-air pressure, the manifold temperature, the manifold humidity, and the manifold pressure; an opening area learning value calculating that calculates a humidity detecting recirculation flow rate which is a flow rate of the recirculation exhaust gas based on the humidity detecting EGR rate and the intake air flow rate, calculating a humidity detecting opening area, which is an opening area of the EGR valve which realizes the humidity detecting recirculation flow rate, and calculates a learning value of the opening area of the EGR valve based on the humidity detecting opening area; and a recirculation exhaust gas calculating for control that calculates a learned opening area of the EGR valve corresponding to the present opening degree of the EGR valve using the learning value of the opening area, and calculates a flow rate of the recirculation exhaust gas for control used for controlling the internal combustion engine based on the learned opening area. 